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ABSTRACT 


In  this  thesis  a  control  systems  analysis  package  is  developed  using  parameter  plane 
methods.  It  is  an  interactive,  user-friendly  computer  aid  to  plot  families  of  performance 
index  cost  curves.  By  superimposing  the  cost  curves  on  the  parameter  plane  curves  the 
designer  is  able  to  choose  values  of  the  parameters  uhich  provide  a  good  compromise 
between  cost  and  dynamic  behavior. 
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I.  INTRODUCTION 


y  When  analyzing  or  designing  a  control  system,  it  is  important  to  study  the  effects 
on  0%'erall  system  performance  of  varying  one  or  more  parameters  (mass,  inertia,  gain, 
resistance,  etc).  It  is  equally  important  to  determine  whether  a  desired  dynamic  behavior 
can  be  achieved  with  any  set  of  values  for  the  parameters  ...  if  not.  redesign  is  indicated. 

The  analysis  and  synthesis  of  linear  feedback  control  systems,  or  the  compensation 
of  same,  can  be  realized  by  parameter  plane  methods.  The  parameter  plane  method, 
which  works  well  for  two  variable  parameters  and  which  may  be  extended  to  three  or 
more  parameters,  is  purely  algebraic,  and  the  resulting  plots  are  valuable  aids  to  analysis. 

Furthermore,  the  performance  index  cost  curves  in  terms  of  integral-square  error 
and  other  cost  functions  can  be  superimposed  on  the  parameter  plane  graphically.- That 
is  veiy  useful  for  analysis  and  for  designing  a  control  system  when  its  characteristic 
equation  has  two  parameters.  Since  the  coetTicients  of  the  characteristic  polynomial  are 
determined  by  the  system  parameters,  it  follows  that  some  relationship  exists  between 
the  value  of  any  parameter  and  the  value  of  the  characteristic  roots.  The  integral-square 
error  as  a  performance  inde.x  is  also  determined  by  the  system  parameters. 

Let  us  e.xamine  the  choice  of  the  integral-square  error  as  a  performance  index.  In 
many  applications  the  performance  of  feedback  control  systems  can  be  considered  either 
satisfactory  if  the  errors  are  below  a  certain  limit,  or  unsatisfactory,  if  the  errors  are 
above  a  certain  lintit.  The  hypothesis  that  the  configuration  is  fi.xed  refers  to  the  form 
of  the  system;  it  is  assumed  that  one  or  more  parameters  can  be  adjusted.  Given  the 
input  and  ideal  output  signals  as  transient  functions  of  time  and  given  the  configuration 
of  the  control  system,  find  the  values  that  the  free  parameters  should  have  in  order  to 
minimize  the  integral-square  error. 

One  approach  to  solving  this  problem  starts  with  the  Fourier  transforms  of  the  input 
and  ideal  output  signals.  The  integral-square  error  is  then  easily  formulated  by  means 
of  Parseval's  theorem.  If  the  transform  of  the  error  is  a  rational  function  of  the  complex 
frequency,  the  integral-square  error  can  be  easily  evaluated  by  means  of  a  table  of  in¬ 
tegral,  which  is  in  .Appendix  A.  When  the  integral-square  error  is  evaluated,  it  appears 
\  as  a  function  of  the  free  parameters. 

'  The  object  of  this  thesis,  then,  is  to  develop  a  user  friendly  interactive  computer 
program  to  plot  farrulies  of  performance  index  cost  curves.  .A  secondary  objective  is  to 
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^  demonstrate  the  use  of  these  curves  in  the  design  of  the  system.  The  design  procedure 
starts  with  the  performance  index  cost  curves  superimposed  on  parameter  plane  curves. 
The  result  of  the  analysis  is  a  choice  of  values  for  the  parameters  such  that  the  system 
performance  is  the  best  available  compromise  between  transient  performance  and  cost. 
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II.  BACKGROUND  THEORY  AND  EQUATIONS 
A.  FUNDAMENTAL  CONCEPT 

Classical  control  theory,  historically  speaking,  has  been  applied  to  the  design  prob¬ 
lem  with  the  assistance  of  graphical  presentation  and  trial  and  error  methods.  Such 
methods  have  been  quite  successful  in  the  development  of  good  control  systems,  but  do 
not  answer  the  question  “Is  this  the  best  system  possible.”  The  engineer  and  mathe¬ 
matician  have  always  looked  for  the  “best”  control  system.  Early  methods  used  “stand¬ 
ard  forms.”  The  concept  of  a  “performance  index”  was  invoked,  also  at  an  early  point, 
for,  if  the  chosen  index  could  be  satisfied,  one  could  say  that  the  resulting  system  was 
indeed  the  “best”  as  measured  by  the  index. 


Let  us  take  the  performance  index  that  is  the  integral-square  value  of  the  error  be¬ 
tween  the  ideal  output  and  the  actual  output.  A  simple  feedback  control  system  with 
adjustable  parameters  is  assumed  as  presented  in  Figure  1. 


Figure  I.  Block  diagram  of  feedback  control  system. 
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The  performance  index  in  terms  of  integral-square  error  is  defined  as; 


PI  = 


jE^d, 


(2.1) 


where  PI  ;  performance  index 

E  ;  input  -  actual  output 

and  the  criterion  was  to  design  the  system  such  that  this  index  was  minimized.  For  lin¬ 
ear  systems  this  problem  can  be  solved  by  applying  Parseval's  theorem  [Ref  1:  p.  43]  to 
obtain  explicit  relationships  for  system  parameters.  One  should  note  that  a  system  de¬ 
signed  to  minimize  this  performance  index  has  a  step  response  with  considerable  over¬ 
shoot  and  oscillation,  and  this  may  not  be  acceptable  for  some  applications,  but  where 
such  dynamic  characteristics  are  suitable  the  design  is  straightforward.  Other  perform¬ 
ance  indices  which  have  been  found  useful  are 


PI=j|£)/^//  (2.2) 

PI=Jl£Ur  (2.3) 

Both  of  these  are  nonlinear  and  do  not  have  analytic  solutions.  Designs  using  such 
criteria  were  accomplished  using  simulation.  The  solution  of  an  optimal  control  problem 
is  the  result  of  minimizing  the  chosen  cost  function  and  is  the  design  which  best  satisfies 
the  criteria  defined  by  that  cost  function. 

B.  PARSEVAL'S  THEOREM 

Suppose  we  wish  to  know  the  area  beneath  the  product  of  two  functions  xM  and 
.Vj(r)  of  the  infinite  interval  of  the  time  from  minus  infinity  to  plus  infinity.  Let  this  area 
be  denoted  by  the  symbol  I  (which  stands  for  integral),  thus 

r+oo 

1=  x,{i)x,{i)di  (2.4) 
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Suppose,  in  addition,  that  the  time  functions  are  Fourier  transformable  and  have  trans¬ 
forms  .v,(s)  and  .Y,(s),  respectively.  We  wish  to  express  the  integral  I  directly  in  terms 
of  these  transforms  in  order  to  obviate  inverse  transforming  .ri(s)  and  .r,(s)  when  these 
are  given  rather  than  the  time  functions  themselves.  To  do  this  we  observe  that 


jC2(s)  ds 


Substitution  of  this  value  of  .r2(r)  into  Equation  (2.4)  yields 


(2.5) 


p+oo  I  C+Joo 

*=J  -^iW  -2~\  e".t2(s)i/s 


dt 


(2.6) 


Let  us  now  interchange  the  order  of  the  integrations  so  that  we  integrate  with  respect 
to  time  first 


I  r  +j'X  r  +CO 

1=^1-  e^'x.iOdi  (2.7) 

By  the  direct  Fourier  transform  we  know  that 

f+oo 

.r,(5)=  e~^'x,(/)Jr  (2.8) 

oo 

Thus  the  integral  with  respect  to  time  on  the  right  side  of  Equation  (2.7)  may  be  evalu¬ 
ated  as 


;c,(0rft  =  .Y,(  -s) 


This  permits  us  to  write  Equation  (2.7)  as 


1  = 


—S)  X2(s)  ds 


This  is  the  desired  expression  for  the  integral  I. 


(2.9) 


(2.10) 
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An  important  particular  case  of  the  above  result  occur  when  the  two  time  functions 
involved  are  identical  so  that  the  integral  I  becomes 


We  call  this  integral  the  integral-square  value  of  the  time  function  .v(r). 
By  Equation  (2.10)  we  know  this  can  be  e.xpressed  as 


This  result  is  known  as  Parseval's  theorem.  This  theorem  represents  a  very  convenient 
way  of  expressing  the  integral-square  value  of  a  time  function  in  terms  of  its  transform. 

C.  MINIMIZATION  PROCEDURE 

On  the  basis  of  Parseval's  theorem  we  are  now  in  a  position  to  formulate  the  pro¬ 
cedure  for  minimizing  the  integral-square  error  for  transient  input  signals  by  adjustment 
of  the  free  parameters. 

This  procedure  involves  four  steps.  The  first  step  is  to  express  the  Fourier  transform 
of  the  error  as  a  function  of  the  complex  frequency  s.  This  function  will  involve  the  free 
parameters  of  the  system  as  unknown  coelTicients.  The  second  step  is  to  express  the 
integral-square  error  I,  in  terms  of  the  error  transform  E(s)  by  Parseval's  theorem,  i.e., 

=  (2.13) 

At  this  stage,  provided  E(s)  is  a  rational  function,  the  integral-square  error  will  appear 
in  the  form 


j  f +y«>  c(  — s)  c(s) 
*  ~  2;ry  d{  -s)  d{s) 


(2.14) 


where  c(s)  and  d{s)  are  polynomials  in  s.  The  third  step  is  to  evaluate  the  integral. 
Fortunately,  definite  integrals  of  this  form  have  been  evaluated  in  terms  of  the  coeffi¬ 
cients  appearing  in  the  polynomials.  A  table  of  these  integrals  along  with  an  explanation 
of  how  they  were  evaluated  is  given  in  Appendix  A.  Upon  carrying  out  the  integration 
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we  have  the  integral-square  error  expressed  as  a  function  of  the  free  parameters 
p,  through  p*.  Symbolically  the  integral-square  error  at  this  stage  is  expressed  as 

^eiPx.Pi,  •••./’fe) 


The  fourth  step  of  the  minimization  procedure  is  to  adjust  the  values  of  the  free  pa¬ 
rameters  in  such  a  way  as  to  minimize  the  integral-square  error.  Formally  this  can  be 
done  in  the  standard  manner  of  equating  the  partial  derivatives  of  the  integral-square 
error  to  zero  and  solving  the  resultant  set  of  equations  for  the  values  of  the  parameters. 
In  symbols  this  means  solving  the  set  of  k  equations  of  the  form 


Unfortunately  in  many  practical  problems  this  formal  procedure  for  finding  the  param¬ 
eter  values  that  minimize  the  integral-square  error  leads  to  sets  of  non-linear  equations 
in  the  parameters  which  are  so  complex  that  straightforward  methods  of  solution  do  not 
exist. 

If  the  equations  for  the  parameters  obtained  by  setting  the  partial  derivatives  of  the 
integral-square  error  with  respect  to  the  parameters  equal  to  zero  are  too  complex  to 
be  solved  by  analytical  methods,  we  may  resort  to  trial-and-crror  or  successive  approxi¬ 
mation  techniques.  However,  if  such  methods  are  necessary  to  solve  the  equations  for 
the  parameters,  it  is  probably  best  to  go  back  to  the  expression  for  the  integral- square 
error  itself  and  plot  this  as  a  function  of  one  of  the  parameters,  holding  the  other  pa¬ 
rameters  fixed.  By  using  a  sufficient  number  of  fixed  values  for  these  other  parameters 
and  making  a  sufficient  number  of  plots,  values  can  be  found  for  all  the  parameters  that 
define  the  minimum  value  of  the  integral-square  error  closely  enough  for  engineering 
purposes.  By  working  directly  with  the  integral- square  error  expression,  we  avoid  taking 
the  partial  derivative  of  this  expression  with  respect  to  each  of  the  free  parameters. 
Taking  these  derivatives  may  be  a  very  difficult  and  time-consuming  chore. 
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III.  PERFORMANCE  INDEX  CURVE  EVALUATION 

A.  BACKGROUND 

1.  Overview 

The  performance  index  program  is  an  interactive  personal  computer  program 
whenever  possible.  Before  using  this  program,  we  have  to  generate  a  datafile  containing 
a  ,  fi  and  PI.  This  file  will  enable  us  to  plot  families  of  performance  index  cost  curves. 
There  currently  exists,  at  Naval  Postgraduate  School,  a  DSL(Dynamic  Simulation  Lan¬ 
guage)  [Ref  2]  program  residing  on  the  school's  IBM  370  mainframe.  By  using  the  DSL 
program,  we  can  easily  obtain  the  above  mentioned  data  file,  and  then,  down  load  this 
file  from  the  mainframe  computer  to  floppy  diskette.  This  data,  now  on  floppy  diskette, 
is  to  be  used,  along  with  the  computer  program  developed  in  this  research  paper,  to  plot 
performance  index  curves.  With  this  program,  the  user  can  examine  particular  values  of 
the  performance  index  curves.  Options,  such  as  printer  type  and  size  of  output,  can  be 
selected  by  the  user. 

2.  Soft^tare  /  Hardware 

The  program  language  used  in  coding  this  program  is  Microsoft  FORTRAN?? 
V3.20.  The  Plotworks  PLOT88  graphics  library  [Ref  3]  is  used  to  generate  output  plots. 
The  plotting  of  performance  index  curves  conducted  for  this  thesis  were  performed  on 
an  IB.M-AT  or  IB.M  compatible  80286  machine.  Due  to  the  universality  of  FORTR.AN 
coding,  this  program  could  be  implemented  on  any  machine  capable  of  being  pro¬ 
grammed  in  FORTR.AN.  The  source  code  of  generating  curves  of  constant  perform¬ 
ance  index  is  listed  in  Appendix  B. 

Before  running  the  program,  ANSI. SYS  should  be  incorporated  in  the  personal 
computer's  CONFIG.SYS  file  prior  to  running  the  performance  index  program. 
Through  the  ANSI. SYS  device  driver,  system  calls  to  clear  the  screen  and  position  the 
cursor  are  enabled.  If  not,  screen  readability  suffers  although  the  program  is  still  fully 
functional.  .More  will  be  said  on  this  in  the  ANSI  module  section  of  this  chapter. 
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3.  T«o  Dimensional  Contouring 


1  he  control  grid  is  constructed  as  follows  Figure  2. 


Figure  2.  The  contour  grid  t»u-diiiiensional  real  array  z  from  Ref.3. 
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where  •  =  dimensional  size  of  real  array  z. 

(nx.ny)  =  maximum  index  of  area  used. 

(xlow.ylow)  =  lower  left  corner  of  the  grid  area  in  the  world  units, 
(xhigh.yhigh)  =  upper  right  corner  of  the  grid  area  in  the  world  units, 
(nxsize.nysize)  =  grid  interval  in  world  units, 
dx  =  (xhigh-xlow)  (nx-1) 
dy  =  (yhigh-ylow),  (ny-l) 


Note  that  only  the  control  points  in  the  range  (i=  1, . ,nx  ;  j  =  1, . ny)  are  used 

in  the  array  z.  By  convention,  the  first  index  is  referred  to  as  the  row  index  and  the  sec¬ 
ond  index  is  referred  to  as  the  column  index. 

The  following  sequence  produces  a  two-dimensional  contour  map: 

•  A  grid  constructed  over  the  surface  Figure  2. 

.At  each  control  point  of  the  grid,  the  z  ordinate  is  determined.  Using  Laplace  or 
Spline  interpolation  or  a  combination  of  the  two.  curves  of  constant  z  are  con¬ 
structed  by  PLOTS8  graphics  librar\'.  To  do  that,  we  can  use  the  subroutine  named 
ZGRID  which  is  the  PLOT88  librarx',  this  subroutine  interpolates  a  set  of  arbi¬ 
trarily  spaced  data  points  xp{k)  ,yp{k),  k=  npis  onto  a  rectangular  grid  z{ij)  of 
arbitrarx'  size. 

•  The  grid  control  points  may  be  smoothed.  The  ZSMTH  subroutine  which  is  the 
PLOT88  library  smooths  the  gridded  control  points  by  applying  Laplacian 
smoothing. 

•  Potions  of  the  contour  may  be  removed  to  make  room  for  annotation. 

•  Tlie  contours  are  drawn.  Each  set  of  4  grid  points  can  be  subdivided  further  into 
4  or  16  subgrids  for  more  precise  contouring.  In  this  program,  the  set  of  4  grid 
points  is  subdivided  into  16  subgrids.  Variable  contour  line  types  and  widths  can 
be  drawn.  To  make  drawing  more  precisely,  we  can  use  the  subroutine  named 
ZeSEG  which  is  the  PLOTSS  librarx-.  This  subroutine  converts  the  grid  of  control 
points  into  a  set  of  contours  and  draws  the  contours.  In  addition,  the  capability 
of  dividing  each  grid  square  into  1,  4,  16  subsections  using  cubic  polynomial  in¬ 
terpolation.  This  step  may  be  needed  if  contours  are  close  together  or  if  the  line 
quality  of  contour  need  improvement. 

The  size  of  the  contour  grid  is  defined  in  the  main  program.  For  most  applica¬ 
tions.  a  size  of  50  x  50  control  points  is  sufficient.  The  maximum  number  of  control 
points  is  limited  by  the  available  space  in  the  data  segment  of  the  program.  A  maximum 
of  the  51  control  levels  can  be  drawn. 
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B.  THE  PROGRAM 


The  performance  index  program  given  in  Appendix  B  itself  consists  of  main  program 
and  six  subroutines  which  make  up  tlie  performance  index  package.  Figure  3  is  a  flow 
chart  illustrating  the  organization. 


Figure  3.  Performance  Inde.x  Program  Flow  Chart. 


When  the  program  is  initiated,  the  following  questions  are  presented; 
1.  IIow  many  constant  PI  curves  ^  =  =  =  > 

•  The  user  has  to  enter  an  integer  value  which  is  must  be  less  than  51. 


2.  Unter  PI  values,  Ivvgt,  Idig  #  =  =  =  > 
•  Table  1  give  us  proper  values. 


Table  1.  PARAMETER  TABLE 


|QQQQQ|||[ 

type 

value 

Description 

PI 

real 

Performance  Indc.x 

hvgt 

integer 

1  to  27 

the  line  type  for  each  PI  curve  the  PI  curve 
is  drawn  with  a  solid  line.  See  PLOT88 
Manual.lRef  3:  p.  32J 

Idit 

integer 

-1  to  -27 

the  PI  curves  is  drawn  with  a  dashed  line. 

Idit 

integer 

>  0 

the  annotation  control  information  for  each 

PI.  number  of  digits  after  decimal  point. 

Idit 

integer 

-1 

omit  decimal  point. 

Idit 

integer 

<  -2 

omit  contour  annotation  for  this  level. 

3.  Enter  number  of  data  file  points  =  =  =  > 

•  The  user  has  to  enter  total  number  of  data  points. 

1.  Subroutine  SIZE 

The  subroutine  SIZE  Figure  4  is  displayed.  This  subroutine  permits  adjustment 
of  plot  size  to  both  the  screen  and  the  printer.  This  option  was  utilized  to  scale  plot 
outputs  for  inclusion  in  this  thesis. 


This  is  the  routine  to  Adjust  the  plot  size  *** 


Output  plots  are  currently  sized  to  fill  the  console  screen. 

Would  you  like  to  adjust  the  plot  size? 

Enter  "y"  or  "n"  ==>  y 

A  default  plot  factor  size  of  0.9  is  used  to  fill  the  console  screen. 
This  plot  size  will  be  halved  by  entering  a  value  of  0.45. 

It  will  be  doubled  with  a  factor  entry  of  1.8. 

Enterfactor  size  ==>  0.7 
Figure  4.  The  plot  size  subroutine. 
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2.  Subroutine  LDFILE 

Ihe  subroutine  LDFILE  (Figure  .5)  is  displayed  on  the  screen. 


***  This  is  the  routine  to  LOAD  data  from  a  file  *** 


File  name  should  not  exceed  8  characters 
File  extension  should  not  exceed  3  characters 

What  is  the  file  name  (fn)  and  extension  (ext)? 

Enter  in  form  fn.ext  (e.g.  PARAPLAN . INP )  ===>  pil.dat 


Figure  5.  Tlie  data  load  subroutine. 

Subroutine  LDFILE  permits  program  input  of  a  pre-existing  data  file.  The  user 
is  queried  as  to  file  name  (not  to  exceed  eight  characters)  and  file  extension  (three  char¬ 
acters  or  less).  The  code  checks  for  existence  of  the  file  within  the  working  subdirectory. 
If  it  dees  not  exist,  a  message  is  returned  stating  such  and  offering  the  option  of  entering 
another  file  name  again.  Entry  of  file  extension  is  optional.  In  fact,  a  file  can  be  labeled 
with  as  little  as  one  character  or  number. 

We  note  that  the  data  should  be  formatted  in  order  to  read  the  data  file  suc¬ 
cessfully  by  using  the  LDFILE.  A  data  file  could  be  created  aligning  in  the  format 
identified  by  the  subroutine  LDFILE  source  code  listed  in  Appendix  B. 

3.  Subroutine  MONPRT 

The  largest  of  the  subroutines  is  subroutine  MONPRT.  A  number  of  smaller 
subroutines  are  associated  with  it.  Subroutine  .MO.NPRT  allows  selection  of  the  graphics 
output  device,  be  it  the  monitor  or  a  wide  variety  of  printers.  Graphics  output  defaults 
to  the  monitor  when  the  program  is  first  executed.  Graphics  output  can  always  be 
dumped  from  the  screen  to  a  printer,  but  there  may  be  times  when  direct  output  to  a 
printer  is  desired.  All  NTS  printers  are  included  in  the  printer  output  options.  In  addi¬ 
tion,  one  option  permits  the  user  to  directly  enter  the  values  for  lOPORT  and  MODEL, 
as  outlined  in  the  PLOT88  manual.  (Ref.  3;  p.  17] 
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Tile  first  menu  displayed  when  subroutine  MONPRT  is  called  is  the 


printer.'output  menu  Figure  6. 


;  PRINTER/OUTPUT  MENU  ! 

m 

PRINTER  NO.  I 

PRINTER 

1 

t 

1 

1 

1  I 

Epson  FX-80,  All 

1 

1 

1 

1 

2  ! 

Epson  FX-lOO,  All 

1 

I 

« 

1 

3  i 

Epson  MX-lOO,  All 

1 

1 

1 

1 

4  ! 

Epson  RX-80,  All 

1 

1 

1 

1 

5  1 

Epson  MX-80  A  IBM  Printer 

1 

1 

1 

1 

6  1 

HP  7470A  Graphics  Plotter 

1 

• 

1 

1 

7  1 

HP  7475A  Graphics  Plotter 

1 

1 

1 

1 

8  ; 

HP  758xB  Series  Plotters 

1 

1 

1 

1 

9  1 

HP  2686A  Laser  Jet 

1 

1 

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

1 

1 

10  ! 

Graphics  Monitor  (default) 

1 

1 

1 

1 

11  1 

HARDWARE  Interface  Menu 

1 

1 

1 

12  I 

Input  PLOT88  Values  for  lOPORT  and  MODEL 

1 

1 

1 

1 

99  ! 

RETURN  to  LOAD  Data 

1 

1 

Enter  integer 

number  for  selection  ===> 

Figure  6.  The  printer/uutput  menu. 


Selection  of  a  printer  in  this  menu  automatically  directs  graphics  output  to  the  most 
commonly  associated  output  port,  either  parallel  or  serial,  for  that  particular  device. 
However,  usually  there  are  multiple  parallel  or  serial  ports  to  which  a  printer  can  be  at¬ 
tached.  The  default  printer  output  setting  may  not  direct  plots  to  the  appropriate  port. 
To  provide  ma.ximum  output  fle.xibility,  the  printer/output  menu  allows  selection  of  a 
specific  output  port  through  access  to  another  menu. 
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HARDWARE  INTERFACE  MENU 
PORT 


SELECT  NO. 


!  1  1 

!  2  ! 

I  O  I 

I  I 

!  4  1  COMl 

I  5  I  COM2 

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 
!  99  1  RETURN  to  LOAD  Data  ! 


LPTl  Printer  Port 
LPT2  Printer  Port 
LPT3  Printer  Port 
Serial  Port 
Serial  Port 


Enter  integer  number  for  selection  ===> 


Figure  7. 


The  Hardware  Interface  Menu. 


The  output  port  selection  option  is  called  the  HARDWARE  INTERFACE 
MENU  Figure  7.  Ihis  menu  is  part  of  a  subordinate  MONPRl'  subroutine  titled 
PORT.  The  option  allow  selection  ofone  of  three  parallel  ports  (LPTl  -  LPT3)  or  either 
of  two  serial  ports  (COMl  -  COM2). 

If  a  serial  port  is  manually  selected  in  the  HARDWARE  INTERFACE  MENU, 
an  associated  data  transfer  (baud)  rate  must  also  be  assigned.  Immediately  after  se¬ 
lection  of  serial  port  options  4  or  5,  the  BAUD  (data  transfer)  R/\TE  MENU  is  dis¬ 
played  Figure  8.  'fransfer  of  graphics  data  over  a  single  line  is  usually  time  consuming, 
therefore  the  highest  printer  capable  transfer  rate  available  (9600  baud)  is  normally  the 
best  selection.  A  slower  transfer  rate  would  only  be  indicated  if  improper  graphics  out¬ 
put  is  generated. 


t 

1 

BAUD  (data  transfer)  RATE  MENU 

;  SELECT  NO. 

:  BAUD  RATE  1 

1  1 

:  300  1 

!  2 

1  1200  ; 

!  3 

;  4800  ; 

!  4 

;  9600  ; 

Enter  integer  number  for  selection  ===> 


Figure  8.  B.iud  (d.at.i  transfer)  Rate  Menu. 
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Bit  by  bit  transfer  of  data  can  also  have  a  check  sum  (parity)  associated  with  it 
to  provide  an  interval  check  of  correct  data  transfer.  Following  selection  of  a  data 
transfer  rate,  the  PARITY  MFNU  Figure  9  is  displayed.  The  user  has  the  option  of 
choosing  odd, even  or  no  parity. 


1  PARITY  MENU 

1  SELECT  NO.  1 

PARITY  1 

!  1  1 

NO  Parity  1 

:  2  ; 

EVEN  Parity  I 

!  3  I 

ODD  Parity  j 

Enter  integer  number  for  selection  ===> 


Figure  9.  Parity  Menu. 

The  PLOTS8  Software  Library  Reference  Manual  (Ref.  3:  p.  18]  contains  a 
number  of  tables  which  provide  conunon  output  settings  for  many  more  printers  than 
can  be  selected  in  the  PRIN'I  ER.  OUTPUT  menu. 

C.  ANSI  MODULE 

Hie  ANSI  module  contains  two  subroutines.  Subroutine  CLR  clears  the  monitor 
screen.  It  is  used  primarily  as  a  precursor  to  displaying  a  menu  or  positioning  queries 
or  informational  tc.xt.  Subroutine  CUP  positions  the  cursor  inunediatcly  after  a  request 
lor  response.  This  subroutine  serves  not  only  to  give  a  pleasing  screen  display  but  also 
provide  the  user  with  an  inunediate  reference  as  to  where  he  is  the  program  and  if  input 
is  required  of  the  user.  Both  routines  can  be  found  in  the  PLO  I  SS  Reference  Manual 
[Ref  3].  The  two  subroutines  which  comprise  the  ANSI  module  make  use  of  ANSI  es¬ 
cape  sequences.  These  sequences  arc  a  series  of  characters  that  can  be  used  to  define 
lunctions  to  MS-DOS.  In  addition  to  clearing  the  screen  and  positioning,  moving,  and 
saving  the  cursor's  position,  ANSI  escape  sequences  can  control  screen  graphics  and 
reassign  key  definitions. 
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In  order  to  utilize  the  escape  sequences,  the  microprocessor  must  first  be  able  to 
interpret  the  commands.  It  does  this  through  a  systems  file  entitled  .ANSI. SYS. 
ANSI. SYS  must  be  defined  as  a  device  in  the  CONFIG.SYS  file  to  enable  proper  inter¬ 
pretation  of  ANSI  escape  sequences  (e.g.,  device=  \  for  \  ansi.sys  where  \  for  \  is  the 
path  to  ansi.sys).  If  ANSI.SYS  is  not  resident  in  the  CONFIG.SYS  file,  the  perform¬ 
ance  index  program  will  not  clear  the  screen  or  position  the  cursor  when  so  directed. 
Instead,  the  indicated  escape  code  sequence  will  be  printed  on  the  monitor  whenever  a 
call  is  made  to  CLR  or  CUP.  Although  program  readability  is  degraded,  the  perform¬ 
ance  index  program  will  still  function  correctly. 
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IV.  PARAMETER  PLANE  DESIGN  -GRAPHICAL  METHOD 

A.  OVERVIEW 

The  algebraic  solutions  have  the  disadvantage  that  a  fixed  value  of  zeta  and  omega 
must  first  be  chosen  to  compute  the  alpha  and  beta  terms.  In  some  instances  it  is  pos¬ 
sible  to  modify  the  remainder  polynomial  so  as  to  ensure  that  the  specified  roots  are 
dominant.  However,  it  is  not  always  possible  to  guarantee  that  roots  placed  at  a  speci¬ 
fied  location  can  be  made  dominant.  Thus,  an  exhaustive  trial-and-error  procedure  may 
be  required  to  achieve  the  best  values  for  the  various  parameters.  To  avoid  this  trial- 
and-error  analytical  approach  one  can  employ  the  graphical  solution. 

We  have  two  programs  available.  One  permits  us  to  get  constant  zeta  curves  and 
omega  curves. The  second  program  gives  us  a  families  of  constant  performance  index 
curves  in  the  parameter  plane.  We  can  incorporate  the  results  of  two  programs  by 
superimposing  the  plotted  performance  index  cost  curves  over  the  parameter  plane 
curves  on  the  same  sheet. 

To  do  that,  we  perform  a  preliminaiy  procedure.  This  procedure  involves  four  steps. 
The  first  step  is  to  plot  the  parameter  plane  curves  for  the  specific  problem.  This  curve 
permits  us  to  choose  ranges  of  two  parameters  that  are  to  be  used  in  the  analysis  and 
design.  The  second  step  is  to  obtain  a  data  file  containing  a  .  /?  and  performance  index 
values  by  using  the  D.S.L  program.  We  then  copy  this  data  file  from  the  mainframe  to 
a  floppy  diskette.  The  third  step  is  to  plot  the  families  of  performance  index  curves  in 
the  same  range  of  parameter  values  on  the  two  parameter  plane.  We  then  have  two 
kinds  of  parameter  plane  curves.  The  fourth  step  of  the  analysis  for  a  given  specific 
problem  is  to  duplicate  these  curves  on  the  separate  transparent  sheets  and  then  put 
them  together  on  the  same  sheet.  Recapitulation  of  the  graphical  analysis  path  is  shown 
in  Figure  10. 


18 


Figure  iO.  Graphical  analysis  path. 


B.  GRAPHICAL  SOLUTIONS 

In  tills  section,  three  diircrcnt  kinds  of  control  problems  will  be  demonstrated.  One 
of  these  is  a  cascade  compensator  system.  The  other  one  is  a  velocity  and  acceleration 
feedback  compensator  system.  The  above  two  problems  will  be  evaluated  with  three 
kinds  of  performance  inde.x  constraints,  i.e.,  Pl  =  fE?dl  ,  ri  =  f\E\i(il  and 
PI  =  f\E\  lit.  fhe  third  problem  is  a  steering  control  for  optimized  autopilots.  The 
controller  parameters  are  adjusted  to  minimize  the  cost  function. 

I  Jq 

where  <5  =  rudder  angle 
i//,=  yaw  error 
/i.  =  a  weighting  factor 
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1.  Example  4-1  (cascade  compensator) 

Figure  1 1  gives  the  block  diagram  of  a  high  gain  third  order  system  which  is  to 
be  damped  with  a  cascade  filter.  {Ref.  4;  p.  319j 


Figure  1 1.  The  third  order  cascade  compensated  system. 


Let 


P 

“  Z 
f}  =  P 

The  cliaracteristic  equation  for  this  system  is  then; 

5^  +  (ll  -l-/i)SV(10+  ll/?)5^-l-(10/i-M00a)5+  100/?  =  0  (4.1) 

Separating  the  coefllcient  terms. 

•  Constant  coeflkicnt:  S'*  +  I  IS’  -I-  lOS’ 

•  ALPHA  coefficient:  1005 

•  BETA  coefTicient:  5’  -H  I  IS’  -I-  lOS  -blOO 

Enter  each  coeflkicnt  array  value  and  desired  zeta  and  omega  values  into  the  computer 
program.  Figure  12  provides  a  look  at  the  combined  zeta  and  omega  graph  spanning  a 
frequency  range  from  1.0  to  30,  (  =  1.0 , 0.5 , 0.0  and  cu„  =  5.0 
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By  inspection  of  tlie  two  parameter  range,  the  alpha  values  vary  from  0.0  to  140. 
and  beta  values  also  vary  from  0.0  to  100.  Now,  within  these  ranges,  simulation  in 
D.S.L  is  used  to  evaluate  the  cost  at  a  large  number  of  points  on  the  parameter  plane. 
The  D.S.L  simulation  program  and  its  resulting  data  file  containing  a  ,  /?  and  three  dif¬ 
ferent  criterion  cost  values  are  given  in  Appendix  C. 

Once  we  have  obtained  a  data  file,  we  download  this  file  from  the  mainframe 
computer  to  a  lloppy  diskette.  To  run  the  performance  index  program,  type  in  the  exe¬ 
cutable  file  name,  COSTFUN.  The  IBM  microcomputers  in  the  NFS  Controls  Lab  are 
set  up  under  the  file  management  system  IDIR.  To  execute  the  program  under  the  this 
system,  enter  the  appropriate  subdirectory  by  using  the  up, 'down  arrow  keys  to  position 
the  cursor  adjacent  to  the  subdirectory  name,  then  depress  the  ENTER  key.  Use  this 
same  procedure  to  select  COSTEUN.EXE  and  push  ENTER  to  run  the  program.  The 
preparatio*.  step  for  the  plotting  desired  curves  will  appear. 

Enter  proper  values  after  each  prompt  as  follows: 

1. How  many  constant  PI  curves  #  ===>6 

2.  Enter  PI  values  ,lwgt,ldig  #  l==s>0 . 27 , I , 2 

2. Enter  PI  values  .Iwgt.ldig  #  2===>0.2,1,1 

2. Enter  PI  values  ,lwgt,ldig  #  3==a>0 . 15 , 1 , 2 

2. Enter  PI  values  ,lwgt,ldig  #  4==  =  >0 . 12 , 1 , 2 

2. Enter  PI  values  ,lwgt,ldig  #  5===>0 . 1 1 , 1 , 2 

2.  Enter  PI  values  ,lwgt,ldig  #  6s=s>0. 1,1,1 

3.  Enter  number  of  data  file  points  sss>86 


Next,  adjust  the  plot  size. 


***  This  is  the  routine  to  Adjust  the  plot  size 


Output  plots  are  currently  sized  to  fill  the  console  screen. 

Would  you  like  to  adjust  the  plot  size? 

Enter  "y"  or  "n"  ==>  y 


A  default  plot  factor  size  of  0.9  is  used  to  fill  the  console  screen. 
This  plot  size  will  be  halved  by  entering  a  value  of  0.45. 

It  will  be  doubled  with  a  factor  entry  of  1.8. 

Enter  factor  size  ==>  0.7 


The  last  step  is  to  enter  a  data  file  name. 
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♦**  This  is  the  routine  to  LOAD  data  from  a  file  ♦** 


File  name  should  not  exceed  8  characters 
File  extension  should  not  exceed  3  characters 


What  is  the  file  name  (fn)  and  extension  (ext)? 

Enter  in  form  fn.ext  (e.g.  PARAPLAN. INP )  ===>  pil.dat 

After  displaying  the  graph  on  the  monitor,  the  program  prompts  you  for  a 
hardcopy  of  the  graph. 

WOULD  YOU  LIKE  HARDCOPY  OF  THESE  PLOTS?  (Y  or  N) 

y 


PRINTER/OUTPUT  MENU 


PRINTER  NO. 


1 
2 

3 

4 

5 

6 

7 

8 
9 

►♦++++++++++++++ 


PRINTER 


Epson  FX-80,  All 
Epson  FX-lOO,  All  ' 

Epson  MX-lOO,  All 
Epson  RX-80,  All 
Epson  MX-80  &  IBM  Printer 
HP  7470A  Graphics  Plotter 
HP  7475A  Graphics  Plotter 
HP  758xB  Series  Plotters 
HP  2686A  Laser  Jet 
++++++++++++++++++++++++++++++++++++++++++++++++ 


10 

11 

12 

99 


Graphics  Monitor  (default) 

HARDWARE  Interface  Menu 

Input  PLOT88  Values  for  lOPORT  and  MODEL 
RETURN  to  LOAD  Data 


Enter  integer  number  for  selection  ===>  5 

Reenter  the  data  file  name  for  the  hard  copy.  Figure  13  is  a  performance  index 
curve  for  constraint  fE?dt.  In  like  manner,  the  other  two  nonlinear  criterion  of  per¬ 
formance  index  in  terms  of  f\E\  t  dt  and  f\E\  di  can  also  be  plotted  by  running  the 
same  program.  Figures  14  and  15  are  shown  for  desired  constraint  J\E\  t  dt  and 
j\E\  di,  respectively. 
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Let  us  consider  the  two  kinds  of  parameter  plane  curves,  Figure  12  and 
Figure  13.  By  superimposing  one  over  the  other,  the  curves  can  produce  the  very  in¬ 
teresting  graph  of  Figure  16.  In  like  manner.  Figure  12  curves  superimposed  on  the 
other  cost  contour  are  shown  on  Figure  17  and  Figure  18.  These  graphs  contain  three 
kinds  of  constant  families  curves;  C  ,  and  cost  curves.  At  this  point,  from  Figure  16, 
we  could  select  arbitrar>’  operating  points.  The  selected  ALPH.\  BETA  pairs  are: 

A  point:  a  =  25.  /?  =  SO 

B  point;  a  =  19.5,  /?  =  60 

C  point;  a  =  9.0,  /?  =  23 

D  point:  a  =  100,  /?  =  70 

Now,  go  back  to  the  characteristic  equation  Eq.  (4.1)  and  substitute  in  this  equation 
each  .ALPHA  BET.A  pair  and  then,  simulate  the  time  response  for  each  one.  The  tab¬ 
ulated  simulation  results  are  shown  in  Figure  19  and  Table  2  on  page  32.  From  the 
simulation  results,  we  find  that  operating  point  D  (a  =  25,  /i  =  SO)  had  the  fastest  setting 
time,  the  lowest  overshoot  and  the  lowest  cost  inde.x.  The  cost  associated  with  this  point 
are: 


^E^dt  =0.11 

Jl£l  tdi  =0.12 

jl£Ur  =0.23 

And  step  response  is  given  in  Figure  19.  In  this  case,  a  cascade  filter  was  designed  with 
the  zero  equal  to  3.2  (  =  80,  25)  and  the  pole  equal  to  80,  using  a=  P,  Z  and  /?  =  P. 
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Figure  16.  Cuiubine  zeta,  omega  and  cost  curves  for  f  E?  dt. 
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Figure  17.  Combine  zeta,  omega  and  cost  curves  for  f\E\  t  dt. 
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Figure  18.  Combine  zeta,  omega  and  cost  curves  for  f\  fl  dt. 
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Table  2.  SIMULATION  RESULTS. 


op¬ 

erating 

point 

parameters 

I’l  values 

M,, 

t. 

a 

P 

fBdt 

f\E\  tdt 

/l£l  dt 

A 

25 

80 

0.27 

0.43 

0.55 

31 

2.3 

li 

19.5 

60 

0.27 

0.43 

0.53 

31 

2.2 

C 

9 

23 

0.3 

0.53 

0.6 

34 

1.9 

D 

100 

70 

O.ll 

0.12 

0.23 

29 

1.3 

The  exact  values,  however,  will  be  chosen  to  achieve  a  compromise  between 
transient  performance  and  cost,  based  on  given  system  specifications.  It  is  of  interest 
to  compare  the  performance  index  between  the  constraint  f  B  di  and  the  constraints 
f\E\  idi  and  / 1  £l  dt  dii  the  same  operating  points. 


2.  Example  4-2  (feedback  compensator) 


This  example  deals  with  a  plant  incorporating  both  velocity  (tachometer)  and 
acceleration  feedback.  The  block  diagram  of  such  a  system  is  depicted  in  Figure  20.  [Ref. 
4:  p.  331] 


Figure  20.  Velocity  and  acceleration  feedback  compensated  system. 
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The  characteristic  equation  for  this  system  is: 


+  (3  +  S0K2)S^  +  (2  +  50A.T)S  +  50  =  0  (4.2) 


Let 

a  =  a:i 

(i  =  K2 

Using  the  same  procedure  as  the  previous  example.  Figure  21  provides  the  combined 
zeta  and  omega  graph  spanning  a  frequency  range  from  1.0  to  30. 

C  =  0.45 

a,„  =  2.24 

This  graph  shows  the  two  parameters  range,  the  alpha  values  vary  from  0.0  to  1.6  and 
beta  values  var}‘  from  0.0  to  1.0. 

Now.  we  evaluate  the  cost  at  a  large  number  of  points  on  this  range  by  using 
D.S.L.  The  D.S.L  simulation  program  and  its  resultant  data  file  containing  a  ,  /?  and 
three  different  cost  values  are  given  in  Appendix  C.  Once  we  have  the  data  file,  we  can 
produce  the  family  of  performance  index  curves.  These  are  shown  in  Figure  22, 
Figure  23  and  Figure  24,  respectively. 


Figure  21.  Constant  zeta  and  omega  curves  for  example  4.2 
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Figure  23.  The  performance  index  curves  for  / 1  £|  t  dt. 
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Let  us  consider  the  parameter  plane  curves,  Figure  25,  Figure  26  and  Figure  27.  These 
were  produced  by  superimposing  Figure  21  over  each  cost  contour  curve.  Figure  22, 
Figure  23  and  Figure  24. 

Figure  25  consists  of  three  families  of  curves:  zeta,  omega  and  cost  curves.  We 
could  select  arbitrary-  operating  points  given  by; 

A  point:  a  =  0.45,  /?  =  0. 1 8 
B  point:  a  =  0.60,  /?  =  0.40 
C  point:  a  =  0.6,  /?  =  0.067 
D  point:  a  =  0.8,  =  0.22 
E  point:  a  =  0.4,  y?  =  0. 125 

The  operating  point  C  (a  =  0.6,  /?  =  0.067)  permits  us  the  lowest  cost  index.  The  costs 
associated  with  this  points  are; 


\E^dt  =  0.47 
Jl£l  tdt  =0.25 
1 1  £Ur  =  0.68 

However,  settling  time  was  too  slow,  r,  =  0.6  sec.  In  this  system,  operating  point  C 
(a  =  0.6,  /?  =  0.067)  determines  the  parameters  K1  and  K2, 

K\  =  0.6 

K2  =  0.067 
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Figure  26.  Combine  zeta,  omega  and  cost  curves  for  f]E\  t 


Figure  27.  Combine  zeta,  omega  and  cost  curves  for  f\E\  dt. 


RESPONSE  IflLPHfl-0.6.BETfl-0.4)  STEP  RESPONSE 


Table  3.  SIMULATION  RESULTS. 


RB 

parameters 

IM  values 

a 

P 

fpd, 

/I  El  tdi 

/|£l  dt 

A 

0.45 

0.18 

0.6 

1.9 

1.1 

21 

4.1 

B 

0.6 

0.4 

0.76 

2.5 

1.4 

20 

5.9 

C 

0.6 

0.067 

0.47 

0.25 

0.68 

21 

6.0 

D 

0.8 

0.22 

0.62 

0.74 

0.95 

1 

2.1 

E 

0.4 

0.125 

0.6 

2.0 

1.1 

23 

4.5 

Table  3  and  Figure  28  shoAv  simulation  results  for  operating  points  A  through  E.  The 
designer  can  then  choose  the  operating  point  (a,  P)  depending  on  his  desired  specifica¬ 
tions. 

3.  Example  4-3  (optimization  of  controller) 

As  a  final  engineering  example,  consider  the  steering  control  for  a  surface  ship. 
The  block  diagram  of  such  a  system  shown  in  Figure  29. 
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Given  that  the  cost  function  is  defined; 


J  =  yJ  + 


where  6  =  rudder  angle 
ij/,  =  yaw  error 
/  =  a  weighting  factor 

In  this  study,  we  used  values  given  by  R.E  Reid  [Ref  5]  for  the  SL-7  container 
ship.  Reid  uses  a  second  order  Nomoto  model  for  the  SL-7, 


G{s)  = 


s{sT  -F  1) 


and  also  uses  a  controller 


G,(5)  = 


Kl{sT\  +  1) 
(577+1) 


His  results  are  given  in  Table  4. 


Table  4,  REID  RESULTS. 


ship 

plant 

weisht- 

controller 

speed 

Knots 

K 

T 

ing  factor 

K1 

T1 

T2 

16 

0.1084 

90.36 

16.796 

0.4556 

89.51 

10.06 

Let  us  assume  T1  and  T2  are  adjustable  free  parameters. 


r  r  0.1084  0.4556(571  +  1) 

"  5(90.365+  1)’  (5  77+  1) 

Let  a=-^- 

The  characteristic  equation  for  this  system  is: 

90.365^  +  (1  +  90.36/?)S^  +  ()?  +0.049a)5  +  0.049/?  =  0 
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The  parameter  plane  curves  examined  within  the  frequency  range  from  O.Ol  to  0.1 
C  =  0.5,  1.0  window  size;  0  ^  ^  35.,  0  ^  ^  1  Using  the  parameter  plane  program,  we 

get  the  praph  of  Figure  30. 


Figure  30.  Constant  zeta  curves  for  example  4.3 
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In  order  to  evaluate  a  family  of  cost  curves,  we  have  to  obtain  a  data  file  on  the  same 
range  of  parameters,  as  Tigure  30  by  simulation.  The  D.S.L  simulation  program  can 
be  implemented  with  ease  by  rearranging  the  block  diagram  in  Figure  31. 


Figure  31.  Rearrangement  of  block,  diagram. 
The  controller  block  becomes 


(5  +  ir) 

is  +  P) 

The  D.S.L  simulation  program  and  data  file  for  this  example  are  given  in  APPENDIX 
C.  Applying  the  previously  discussed  procedures  for  producing  constant  cost  curves 
provides  the  graph  of  Figure  32. 


46 


Figure  32.  The  cost  curves  for  -jr  dt. 
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From  Figure  32.  we  could  select  the  operating  points  given  by: 


A  point;  a  =  9,  /?  =  0. 1 
B  point:  a  =  20,  /?  =  0.2 
C  point:  a  =  20,  ^  =  0.4 

The  operating  point  A  (a  =  9.0,  j5  =  0.1)  provides  the  minimum  cost. 
Recall  that 


a  = 


IL 

T2 


The  two  parameters  values  become 


P  = 


T2 


n  =  90.0,  T2  =  10.0 


Comparison  of  Reid's  result  and  the  simulation  result  show  that  they  are  very  closely 
matched.  In  addition,  we  simulated  yaw  response  and  rudder  activity  at  the  operating 
points  A  ,  B  and  C  in  Figure  33.  Figure  34  compares  the  yaw  response  to  course 
change  for  the  ship  at  each  operating  point.  Figure  35  compares  the  rudder  activity. 
It  is  noted  that  operating  point  C  is  a  faster  controller  but  has  a  higher  cost  value.  This 
cost  is  a  higher  rate  of  fuel  consumption  used  to  overcome  the  higher  hydrodynamic 
drag  forces. 
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The  yan-  response  vs  time. 


Figure  35.  The  rudder  activity  vs  time. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  graphical  solution  of  the  parameter  plane  technique  provides  the  designer  with 
a  visual  means  of  deducing  how  the  dominant  roots  of  the  characteristic  equation  move 
about  in  the  s-plane  as  two  user-defined  parameters  defining  compensator  attributes  are 
varied.  This  was  illustrated  in  the  example  problems  in  Chapter  4.  The  graphical 
method  is  able  to  analyze  easily,  not  only  the  linear  constraint  jBdt  but  also  the  non¬ 
linear  constraint  j\E\  i  dt ,  J\E\  di  and  any  kind  of  cost  function  using  the  perform¬ 
ance  index  program. 

Constant  cost  value  contours  can  be  plotted  in  much  the  same  presentation  as  the 
parameter  plane  plot.  This  procedure  is  a  useful  design  tool  which  permits  simultaneous 
analysis  of  both  dynamic  response  and  cost. 

As  has  been  demonstrated,  rapid  solution  of  a  wide  range  of  linear  automatic  con¬ 
trol  systems  problems  can  be  achieved  through  use  of  the  performance  index  program 
developed  as  part  of  this  thesis. 

B.  RECOMMENDATIONS 

Further  work  should  be  done  to  combine  individual  programs  so  that  both  param¬ 
eter  plane  and  cost  curves  are  produced  by  one  consolidated  program. 
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APPENDIX  A.  TABULATED  VALUES  OF  THE  INTEGRAL  FORM 


Table  A.l  gives  the  value  of  I„  for  values  of  n  from  1  to  9  where 


. _ 1_  ,  c(5)c(  -5) 

"  2nJ  dis)d{-s) 

and 

c(5)  =  +  •  •  •  +  Co 

<^(s)  =  d„s''  +  •••  +  df) 


(A.l-l) 

(>4.1-2) 

(^.1-3) 


Table  A.l 


‘  2dod^  ' 

_  C\d^  + 

'  2d^d,d,  ' 

2  2  2 
_  Ci^()d\  +  (C)  —  Ic^c-^d^d-^  +  Co^2^3  . 

^  *"  2dQd^{  -d^dj  +  d^d2) 

j  C3(  —dod^  +  d^d^d^)  +  (C2  2c\C2)di^d^dn  +  (C|  2cqC2)c/q(/3^4  +  Co(  ~d^d^  +  d^d-^d^  ^ 

2</oc(i(  —d^dl  —  d\d^-\-  d^d^d^ 

I5  =  [ci/wo  +  (<^3  -  2c2C4)/m,  +  (C2  -  2c, C3  +  2coCi)m2  +  (cf  -  2coC2)m3  +  Cgm^]; 
where 

Wo  =  -^(i/3m,  -i/imj), 
wi]  =  —d^dj  +  did2, 
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m2  ==  -d^d^  d^d^, 

"»3  =  “T  -  ^4"Jl). 

«o 

m^  =  {djm^  —  and 
da 

As  =  dQid^m^  -  d^m^  +  d^mi). 

h  =  [c5"io  +  (cl  -  2c3C5)mi  +  (cj  -  2C2C4  +  2cjC5)m2  +  (cl  -  2c, C3  +  2coC4)m3] 

+  [(‘^?  “  2coC2)m4Com5]; 

where 

mo  =  —  (d4m,  -  4/2/712  +  4/0/713), 

“6 

2  2 

2 

W2  =  ^o^3^5  ^1  ^6  “ 

/7I3  =  4/04/5  +  4/,  jVg  -  4/,</44/5, 

/7I4  =  -J-  (4/2/7I3  -  4/4/7I2  +  l/g/Tl,), 

“0 

/715  =  -j-  (4/2/714  —  4/4/713  +  4/6/712),  and 
“0 

^6  =  4(^1  "*5  -  4'”4  + 

h  =  [<^6"*o  +  (4  -  2c4C6)/7i,  +  (C4  -  2C3C5  +  2c2Ce)/7i2  +  (cj  -  2C2C4  +  2c, C5  -  2coC^)m2] 

+  [(4  -  2c, Cj  +  2CoC4)/7J4  +  (cf  -  2CoC2)/7l5  +  Cgm^J; 


54 


where 


wio  =  -7-  {d^m^  -  +  d^m^), 

“7 

m]  =  —  (d^dtj^  —  dfjd^)  +  {dod^  —  d^d^{dQd^  —  d^d^  +  d^d^  —  d^d^^, 

^2  ~  i^Q^i  ~  ~^o^5  +  (dod^  —  d^d^{d2di  —  d^d^), 

m2  =  —  (df^dy  —  d^d^)  +  {dQd2  —  d^d2){d^d2  —  d^d^), 

m4  =  -j-  (i/jmj  -  d^mj  +  d(,m^), 

“0 

"»5  =  “T  -  ^4"i3  +  d6m2), 

«o 

m^  —  -~  (<72^5  ~  <^4'”4  +  <^6"*3)i  ^nd 
“0 

A7  =  i/oW'”6  -  +  ^5'«4  -  ^7"«3)- 

“  1^7  '''  ~  ^  ^^5  -  2C4C6  +  2C3C7)W2  +  (ci  -  2C3C5  +  2C2C6  -  20,07)^3] 

+  [(C3  -  2C2C4  +  2C1C5  -  2coC6)m4  +  (C2  -  2C1C3  +  2coC4)m5  +  (cf  -  2coC2)m6  +  Com7]; 

where 


Wq  =  —  (fl'gW]  -  +  ^4^3  -  </oW4), 


^”1  ~  (^0^7  ~  ^2^5)1  ~^o^!^7  +  dod2d$  +  ^d\d^)  +  (^^3^/7  —  d\){dQd^  +  d^d\) 

+  did2dg{dQd2  —  d\d2)  —  d^d^^dQd^  —  d^dj) 

4-  (  —<^2^7  +  ^3^6  ~  ^4'^5)(‘^0‘^3  ^1^4)  “  +  l>dQd2d^ 

—  d]d2d2{d2d(,  —  </4</5)  +  2d^^d^d2, 


2 

^2  ~  (^0*^3  ~  ^i^2)(‘^^7  ~  d\dsdi  —  d^df^d^  +  d2d2d2) 
+  {d^d^  —  did2){  —d^d^d^  +  d^d^  —  </i</2^3  +  d\d^) 
—  d^d^d^id^d^  —  +  c/)  d^{dQdj  —  </]<4)» 
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nil  —  d-^id^d^  —  d^d-j)  +  (  —d^d^  +  d(^d-i){d(^d\d^  —  d^d-^  +  d'^d^d^  —  d\d^ 
4"  ^0^7 (  “"^0^5  d-^d^  "h  d^d^  —  'id^d^d'jd^^ 


^4  —  (  —d^d^  +  df/ii){2dQdyd-^  —  d^^d^d^  +  d^d^d^  —  d\d(^) 
+  (  —  d^d^  +  d^d-^Xd^d^d-j  —  d^d^df  +  d\d^  —  d^d^. 


"»5  =  "T  {d2'^4  -  <^4>^3  +  d6m2  -  dg/Tt^), 
“0 


"*6  =  “T  -  4}W4  +  4^3  -  dgin^), 

“o 


m-,  —  ~  {d2m(^  —  d^mg  +  d(,mi  —  dgm^).  and 
“o 


Ag  =  d^id^m-,  -  dgrrtf,  +  d^m^  -  d^m^. 


h  =  [c8"\)  +  (c?  -  2c6C8)m,  +  {cl  -  2csC7  +  2c4Cg)m2  +  (c^  -  2C4C6  +  2C3C7  -  2c2Cg)m2] 

+  [(C4  -  2C3CS  +  2C2C6  -  2c, C7  +  2coCg)m4  +  (cl  -  20^04  +  2c,Cg_2coC(,)mg] 

+ [(^2  -  2C1C3  +  2coC4)m6  +  (Cj  -  2coC2)m7  +  CQmg]; 


where 


^o  =  -7-  (^7"^)  -  d5in2  +  dgmg  -  d^m^). 


—  a,(a,a,o  —  +  OgOg  +  2040^ 


2 

«5 


«5«7  -  «?  ) 


+  «2(  -  02^6  -  «3«7  +  «4«5  +  '2-aA^l)  “  oj, 


^2  ~  ^i(‘23^J9  +  ^4^9  ~  ^5<2g  +  ~  <*?%)  "h  ^2^  ~  ^2^  ~  ^4^7, 


2  2 
"*3  ~  ^i(^3^io  +  «4®io  +  ^7^9  ■"  ^^g)  +  OjC  —  ^2^10  +  ‘27ag)  —  , 


=  a,(a5a,o  +  2a7a,o  -  a8a9)  +  a2(«7«9  "  «4«io)  “  «7. 


Wj  =  -~  {d2m4  ~  d^mg  +  4^2  -  dgm^). 
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mg  =  -j-  {d^m^  -  d^nti  +  d(,m^  -  d^m^), 

«o 

m?  =  ^  (4mg  -  d^m^  +  d(,mi  -  tj'gmj), 

«o 

mg  =  -}-  (i/2m7  -  d^n^  +  df^m^  -  d^m^),  and 

«o 

A,  =  doid^nti  -  t/3m7  +  t/jmg  -  d^ms  +  d^nii), 

where 

—  ^0^3 

^2  ~  ~  ^0^5 

dj  =  <!/3</4  —  </2<^s 

^4  —  did(,  —  d^d-j 
—  d^d(j  —  d^d-i 
Og  =  —  d^di 

a?  =  d^d^  —  dQdg 

ag  =  i/gd’g  —  <^2^9 

a^  —  d^d^  —  d^d^ 

“^10  ~  d-jd^  —  <4^9- 
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APPENDIX  B.  PERFORMANCE  INDEX  PROGRAM. 


This  appendix  contains  a  listing  file  of  the  FORTRAN  source  code  for  the  per¬ 
formance  index  program.  This  program  is  used  to  plot  the  performance  index  curves  for 
desired  values. 

The  program  consists  of  main  program  and  subroutines.  The  subroutine  SIZE  per¬ 
mits  adjustment  of  plot  size  to  both  the  screen  and  the  printer.  The  subroutine  LDFILE 
permits  program  input  of  a  pre-existing  data  file.  Subroutine  MONPRT  allows  selection 
of  the  graphics  output  device,  be  it  the  monitor  or  a  wide  variety  of  printers. 


******  ARRAY  DEFINITION  ****** 


CHARACTER*!  CHAR.CHG 
REAL  Z(51,51) 

REAL  XP(202),YP(202),ZP(200) 
REAL  2PIJ(200) 

INTEGER  KNXT(200),ITEXT(1) 
INTEGER  J.K.SEL 
REAL  ZLEV(50) 

INTEGER  LDIG(50),LWGT(50) 
COMMON/ PLOTTR/  lOPORT, MODEL 


******  VARIABLE  INITIAL  VALUE  ****** 


DATA  XMAX/11.  0/,YMAX/8.  0/ 

DATA  XLPLOT/1. 0/ ,YLPLOT/0.  5/,XHPLOT/8. 0/ ,YHPL0T/5. 5/ 
DATA  CAY/5.  0/,NRNG/3/ 

DATA  HGT/. 12/,NDIV/4/,NARC/8/ 

DATA  NSM/4/ 


******  enter  input  parameters  ****** 

10  CONTINUE 
CALL  CLR 
CALL  CUP(1,0) 

MODEL=99 
IOPORT=99 
13  WRITE(*,297) 

WRITE(*,305) 

READ(*,*)  NLEV 

IF  (NLEV.  LT.  0.  OR.  NLEV.  GT.  50)  GO  TO  13 
DO  22  K=1,NLEV 
WRITE(*,304)  K 

READ(*,*)  ZLEV(K),LWGT(K),LDIG(K) 

22  CONTINUE 

WRITE(*,300) 

READ(*,*)  NPTS 
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WRITE(*,297) 


*  *  *  *  *  *  BEGIN  PROCESSING  *  *  *  *  * 


FACT=0.  9 

CALL  SIZE  (FACT) 

25  CONTINUE 

CALL  LDFILE(NPTS,XP,YP,ZP) 

C.  .  .  INITIALIZE  PLOT88 

CALL  PL0TS(0,I0P0RT, MODEL) 

C. . .  DEFINE  OUTPUT  WINDOW 

CALL  WINDOW(0.  ,0.  ,XMAX,YMAX) 

C. . .  REORIGIN  THE  ENTIRE  DRAWING 
CALL  PL0T(1. 0,0. 5,-3) 

XPMAX=XP( 1) 

YPMAX=YP(1) 

XPMIN=XP(1) 

YPMIN=YP(1I 
DO  3  K=2,NPTS 

IF  (XP(K).GT.  XPMAX)  THEN 
XPMAX=XP(K) 

END  IF 

IF  (YP(K).GT. YPMAX)  THEN 
YPMAX=YP(K) 

END  IF 

IF  (XP(K).  LT. XPMIN)  THEN 
XPMIN=XP(K) 

END  IF 

IF  (YP(K).  LT.  YPMIN)  THEN 
YPMIN=YP(K) 

END  IF 

3  CONTINUE 

DO  4  J=1,NPTS 

XP( J ) =XP( J ) *XHPLOT/XPMAX 
YP ( J ) =YP ( J ) *YHPLOT/ YPMAX 
XHIGH=XHPLOT 
YHIGH=YHPLOT 
XLOW=XPMIN*XHPLOT/XPMAX 
YLOW=YPM I N*YHPLOT/ YPMAX 

4  CONTINUE 

C. . .  ZERO  GRID  ARRAY  Z 
DO  200  1=1,51 
DO  200  J=l,51 
Z(I,J)=0.  0 
200  CONTINUE 

C. . .  GRID  SIZE  PARAMETERS 
NX=8 
NY=6 
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DX=( XPMAX-XPMIN) /FLOAT( NX- 1 ) 

DY=( YPMAX-YPMIN) /FLOAT( NY- 1 ) 

CALL  FACTOR  (FACT) 

CALL  STAXIS  (. 13,. 20,. 15,0. 1,2) 

CALL  AXIS(1.0,0.5,'Alpha',-5,-7.  ,0.  ,XPMIN,DX) 

CALL  AXIS(1.  0,0.  5,'Beta',4,-5.0,90.  ,YPMIN,DY) 

CALL  ZGRID(Z,51,51,NX,NY,XLOW,YLOW,XHIGH,YHIGH,XP,YP,ZP, 

1  NPTS,CAY,NRNG,ZPIJ,KNXT) 

C.  .  POST  THE  INPUT  POINTS. 

CALL  ZPOSTC  XLOW , YLOW , XHI GH , YHIGH , XLPLOT , YLPLOT , XHPLOT , YHPLOT , 

1  XP,YP,ZP,NPTS,2, '(lx,f5. 2)  1389) 

IRTN=590 

C . PUT  A  BOX  AROUND  THE  MAP . 

580  CALL  PLOT(XLPLOT, YLPLOT, 3) 

CALL  PL0T( XHPLOT, YLPLOT, 2) 

CALL  PLOT( XHPLOT , YHPLOT , 2 ) 

CALL  PLOT( XLPLOT, YHPLOT, 2) 

CALL  PLOT( XLPLOT, YLPLOT, 2) 

IF(IRTN.  EQ.  900)  GOTO  900 

590  CALL  PL0T(1.  0,0. 5,-999) 

CALL  PLOTCl. 0,0. 5,-3) 

CALL  FACTOR  (FACT) 

CALL  STAXIS  (. 13,. 20,. 15 ,0. 1,2) 

CALL  AXIS(1. 0,0. 5, 'Alpha' ,-5,-7.  ,0.  ,XPMIN,DX) 

CALL  AXIS(1. 0,0.5,'Beta',4,-5.0,90. ,YPMIN,DY) 

C  .  SMOOTH  THE  DATA  BEFORE  DRAWING  . 

CALL  ZSMTH(Z,51,51,NX,NY,NSM) 

CALL  FACTOR  (FACT) 

CALL  STAXIS  (. 13,. 20,. 15,0. 1,2) 

CALL  AXIS(1.  0,0.  5,  'Alpha\-5,-7.  ,0.  ,XPMIN,DX) 

CALL  AXIS(1.  0,0.5,'Beta',4,-5.0,90.  ,YPMIN,DY) 

CALL  ZCSEG( Z , 5 1 , 5 1 , NX , NY , XLPLOT , YLPLOT , XHPLOT , YHPLOT , 

1  ZLE V , LD IG , LWGT , NLE V , HGT , ND I V , NARC ) 

IRTN=900 
GO  TO  580 

900  CALL  PLOT( 1.0,1.  0,999) 

WRITE(*,*)'  WOULD  YOU  LIKE  HARDCOPY  OF  THESE  PLOTS?  (Y  or  N) ' 
READ(*, 401) CHAR 

IF((CHAR.EQ.  'Y').OR.  (CHAR.  EQ.  'y'))  THEN 
CALL  monpr; 

GOTO  25 
ENDIF 

WRITE(*,*)'  WOULD  YOU  LIKE  TO  DO  ANOTHER  ?  (Y  or  N)' 
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READ(*,401)  CHAR 

IF((CHAR.  EQ.  'Y').OR.  (CHAR.  EQ.  'y'))  GOTO  10 


297  FORMAT(lx) 

300  FORMATC '  3.  Enter  number  of  data  file  points  =>' ,  ) 

304  FORMATC  2.  Enter  PI  values  .Iwgt.ldig  #  ',i2,'=>',  ) 

305  FORMAtC  l.How  many  constant  PI  curves  #  =>' ,  ) 

401  FORMATC A) 

C. . .  end  of  program 
STOP 
END 


SUBROUTINE  CLR 

SUBROUTINE  TO  ERASE  THE  DISPLAY. 
CALL,  CALL  CLR 


SUBROUTINE  CLR 
CHARACTER*!  C1,C2,C3,C4 
INTEGER*2  IC(4) 

EQUIVALENCE  (Cl , IC( 1) ) , (C2, IC(2)) , (C3, IC( 3) ) , ( C4, IC(4) ) 
DATA  IC/ 16#1B , 16#5B , 16#32 , 16#4A/ 

WRITE  ESCAPE  CODE  TO  DISPLAY 
WRITE(*,1)  C1,C2,C3,C4 
1  FORMATC lx, 4a 1) 

RETURN 

END 


SUBROUTINE  CUP 

MOVE  THE  CURSOR  TO  LINE  N,  COLUMN  M. 
CALL,  CALL  CUP(N,M) 

N  -  NEW  CURSOR  LINE 
M  -  NEW  CURSOR  COLUMN 


SUBROUTINE  CUP  (N,M) 

CHARACTER*!  Cl ,C2 ,C5 ,C8 ,LC(5) 

CHARACTER*5  CBUFF 
INTEGER*2  IC(4) 

EQUIVALENCE  (Cl , IC( 1) ) . (C2, IC(2)) ,(C5 , IC( 3) ) , (C8, IC(4) ) , 
*(CBUFF,LC(1)) 

DATA  IC/ 16# IB , 16#5B , 16#3B , 16#66/ 

L=10000+100*N+M 

WRITE(CBUFF,2)L 
2  FORMATC i5) 

WRITEr*,!)  C1,C2,LC(2),LC(3),C5,LC(4),LC(5),C8 
1  FORMATC lx, Sal,  ) 

RETURN 

END 


* 
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C  SUBROUTINE  LDFILE  --  ALLOWS  LOADING  OF  DATA  FROM  AN  EXISTING  FILE  * 


SUBROUTINE  LDFILE  (NPTS.XP.YP.ZP) 

REAL  XP(202),yP(202),ZP(200) 

INTEGER  NPTS.I 
CHARACTER*12  FILEIN.ANS*! 

LOGICAL  EXIST 
1  CALL  CLR 
WRITE(*,100) 

100  F0RMAT(1X,////,10X,’ . ''"'y'Z'll 

*  . ' ,/  lOX,'***  THIS  IS  THE  ROUTINE  TO  LOAD  DATA  FROM  A  FILE  ** 

**' ,/,iox, . . . 

*,//,10X,'FILE  name  SHOULD  NOT  EXCEED  8  CHARACTERS  ./, 

*10X,'FILE  EXTENSION  SHOULD  NOT  EXCEED  3  CHARACTERS  ) 

WRITE!*, 101) 

101  FORMAT! IX, ////,10X,' WHAT  IS  THE  FILE  NAME  !FN)  AND  EXTENSION  !EXT) 

*?  '  ,//,10X, 'ENTER  IN  FORM  FN.  EXT  !E.  G.  PARAPLAN.  INP)  ==>  ) 

CALL  CUP! 18,57) 

READ!*, 102)  FILEIN 

102  FORMAT! A1 2) 

INQUIRE!FILE=FILEIN,EXIST=EXIST) 

IF  ! EXIST)  THEN 

OPEN!7,FILE=FILEIN,STATUS='OLD' ,ACCESS=' SEQUENTIAL' ) 

ELSE 

WRITE!*, 103) 

103  FORMAT!///,! OX, 'THAT  FILE  DOES  NOT  EXIST/, lOX, 

*  'DO  YOU  WANT  TO  TRY  AGAIN?  !Y/N)’) 

READ!*,'!A1)')  ANS  ,  , 

IF!!ANS  .EQ.  'Y')  .OR.  !ANS  . EQ.  y  ))  THEN 

GOTO  1 
ELSE 
GOTO  2 
END  IF 
END  IF 

DO  190  I  =  1,NPTS 

READ! 7, 120)  XP!I),YP!I),ZP!I) 

120  FORMAT! 1X,F10.  2,3X,F10.  2,3X,F10.  3) 

190  CONTINUE 

CLOSE ! 7, STATUS^' KEEP' ) 

2  RETURN 
END 


SUBROUTINE  MONPRT  --  ALLOWS  DISPLAY  ON  MONITOR  OR  PRINTER 


SUBROUTINE  MONPRT 
INTEGER  lOPORT, MODEL, SEL 
COMMON/ PLOTTR/  lOPORT, MODEL 

CALL  CLR 
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80 

81 

87 

82 

83 


WRITE(*,80) 

WRITE(*,81) 

WRITE(*,82) 

WRITE(*,81) 

WRITE(*,83) 

WRITE(*,87) 

WRITE(*,84) 

WRITE(*,81) 

WRITE(*,93) 

FORMAT(lx,///,10x  ’ . 

*  . \/,10x,'  2  lx, 'PRINTER/OUTPUT  MENU',  22x,'  ') 

F0RMAT(10x . 

*  . ') 

FORMATC  lOx f-t  IH  I  +  -t'+f-H +4  K  i-H  \  H  I  ■< -fl  H  1  + 

*4-H-H-H-t-H-’  ) 

FORMATClOx,'  ', 2x, 'PRINTER  NO. ’,2x,’  ' ,20x, ' PRINTER' , 19x, '  ') 
FORMATClOx,'  ',3x,'  1  ' ,2x,'  ' ,4x, 'Epson  FX-80,  All' 


*,26x,' 
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*,/,10x,'  ’ 

,3x.' 

2 

CM 

',4x,' Epson  FX-lOO,  All' 

*,25x,'  ' 
*,/,10x,'  ' 

,3x,' 

3 

’,2x,' 

'  ,4x, 'Epson  MX-lOO,  AH' 

*,25x,'  ' 
*,/,10x,'  ' 

.3x,' 

4 

',2x,' 

'  ,4x,  'Epson  RX-80,  AH' 

*,26x,'  ' 
*,/,10x,'  ' 

,3x,' 

5 

',2x,' 

*,4x, 'Epson  MX-80  &  IBM  Printer' 

*,17x.' 
*,/,10x,'  ' 

,3x,' 

6 

’,2x,' 

’,4x,'HP  7470A  Graphics  Plotter' 

*,17x,'  ' 
*,/,10x,'  ' 

,3x,' 

7 

’,2x,' 

',4x,'HP  7475A  Graphics  Plotter' 

*,17x,'  ' 
*,/,10x,'  ’ 

,3x,' 

8 

',2x,’ 

',4x,'HP  758xB  Series  Plotters' 

*,18x,'  ' 
*,/,10x,'  ' 

,3x,' 

9 

’,2x,' 

’,4x,'HP  2686A  Laser  Jet' 

*,24x,’  ') 
formate lOx 
*'  10 

,’  ',3x 
' ,2x,' 

’',4x. 

'Graphics 

Monitor  (default)' 

*,16x,'  ' 
*,/,10x,'  ' 

,3x,' 

11 

',2x,' 

' ,4x,' HARDWARE  Interface  Menu' 

*,19x,'  ' 
*,/,10x,'  ' 

,3x,' 

12 

',2x,' 

',4x,' Input  PL0T88  Values  for  lOPO 

*RT  and  MODEL' 
*,2x,'  ' 

*,/,10x,'  ',3x,' 

99 

' ,2x,' 

',4x, 'RETURN  to  LOAD  Data',23x, 

93 


*'  '  ) 

FORMATC  lx,///, 15x, 'ENTER  INTEGER  NUMBER  FOR  SELECTION  =>  ') 
CALL  CUP(24,56) 

READC*,*)  SEL 
IF  (SEL  .EQ.  1)  THEN 
lOPORT  =  0 
MODEL  =  5 

ELSE  IF  (SEL  .  EQ.  2)  THEN 
lOPORT  =  0 
MODEL  =  15 

ELSE  IF  (SEL  .  EQ.  3)  THEN 
lOPORT  =  0 
MODEL  =  11 
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ELSE  IF  (SEL  .  EQ.  4)  THEN 
lOPORT  =  0 
MODEL  =  1 

ELSE  IF  (SEL  .  EQ.  5)  THEN 
lOPORT  =  0 
MODEL  =  1 

ELSE  IF  (SEL  .  EQ.  6)  THEN 
lOPORT  =  9600 
MODEL  =  20 

ELSE  IF  (SEL  .  EQ.  7)  THEN 
lOPORT  =  9600 
MODEL  =  30 

ELSE  IF  (SEL  .  EQ.  8)  THEN 
lOPORT  =  9600 
MODEL  =  80 

ELSE  IF  (SEL  . EQ.  9)  THEN 
lOPORT  =  9600 
MODEL  =  60 

ELSE  IF  (SEL  .  EQ.  11)  THEN 
CALL  PORT 

ELSE  IF  (SEL  .  EQ.  12)  THEN 
CALL  CLR 
WRITE(*,100) 

100  F0RMAT(lx,////,5x,' INPUT  VALUE  TO  BE  USED  BY  THE  PL0T88  GRAPHICS 

*  PACKAGE' ,/,5x,' FOR  lOPORT  =>',//) 

CALL  CUP(7,21) 

READ(*,*)  lOPORT 
WRITE(*,110) 

110  FORMAT( lx, ///,5x,' INPUT  VALUE  TO  BE  USED  BY  THE  PL0T88  GRAPHICS 
*PACKAGE' ,/,5x, 'FOR  MODEL  =>',//) 

CALL  CUP( 12,20) 

READ(*,*)  MODEL 

ELSE 

lOPORT  =99 
MODEL  =99 

END  IF 

WRITE(*,98)  lOPORT, MODEL 

98  F0RMAT(//,10x,'PL0T88  WILL  USE  THESE  VALUES  TO  OUTPUT  GRAPHICS:  ' 

*,//, lOx,' lOPORT  =  ' ,14, lOx, 'MODEL  =  ’,14,///) 

PAUSE 

RETURN 

END 


* 

SUBROUTINE  PORT  --  ALLOWS  SELECTION  OF  OUTPUT  PRINTER  PORT  * 


SUBROUTINE  PORT 

INTEGER  lOPORT, MODEL, SEL 

COMMON  PLOTTR/  lOPORT, MODEL 

CALL  CLR 
WRITE(*,79) 

79  FORMAT( lx, ///,15x,' SELECTION  OF  A  PRINTER  IN  THE  PREVIOUS  MENU  AUT 
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*OMATICALLY’ ,/,10x, 'SELECTS  THE  MOST  COMMONLY  ASSOCIATED  OUTPUT  POR 
*T,  BE  IT  PARALLEL  ',/,10x, 

*'(LPT)  OR  SERIAL,  FOR  THAT  PARTICULAR  DEVICE.  IF  THE  PROGRAM  GRAPH 
*ICS' ,/,10x, 

*’ARE  NOT  BEING  OUTPUT  CORRECTLY,  USE  THIS  MENU  TO  PROPERLY  ROUTE  T 
*HE' ,/,10x, 

*' GRAPHICS  DATA  TO  THE  OUTPUT  DEVICE. ',////) 

PAUSE 


80 

81 

87 

82 

83 


CALL  CLR 
WRITE(*,80) 

WRITE(*,81) 

WRITE(*,82) 

WRITE(*,81) 

WRITE(*,83) 

WRITE(*,87) 

WRITE(*,84) 

WRITE(*,81) 

WRITE(*,93) 

FORMATC lx,///,10x  ' . 

*  . \/,10x,’  ’  ,19x,‘ HARDWARE  INTERFACE  MENU' 

*  20x  '  '  '1 

FQRMATdOx,' . 

*  - ') 

FORMAT(  lOx  j  '  ++■(  m  I  H  I  1 4  I  i  I  I  f+4  H- 1  i H  \  I  M  MM-4  I  + 1  H  M  M  I  I  I  IH  H  1  t'4 

FORMATClOx,'  ', 2x, 'SELECT  NO.  ',3x,'  ' , 20x, ' PORT' , 22x, '  ') 
FORMATClOx,'  ',3x,'  1  ’,2x,'  ’, 14x, ' LPTl  Printer  Port' 


84 

93 


*,15x,' 

*,/,10x, 

*,15x,' 

’’M5x,' 

*,/,10x, 

’M6x,' 

*,/,10x, 

*,16x,' 


',3x,'  2  ’,2x,’  ’ ,14x,’ LPT2  Printer  Port' 

'j3x,'  3  ’,2x,'  ',14x,'LPT3  Printer  Port' 

’,3x,'  4  ',2x,'  ’,14x,’C0Ml  Serial  Port' 

',3x,'  5  ',2x,'  ',14x,' COM2  Serial  Port' 


) 


FORMATClOx,'  ',3x, 

*'  99  ',2x,'  ' ,14x,' RETURN  to  LOAD  Data' ,13x,'  ') 

FORMATC lx, ///,15x,' ENTER  INTEGER  NUMBER  FOR  SELECTION  = 
CALL  CUP( 21,56) 

READC*,*)  SEL 
IF  (SEL  .EQ.  1)  THEN 
lOPORT  =  0 

ELSE  IF  (SEL  . EQ.  2)  THEN 
lOPORT  =  2 

ELSE  IF  (SEL  .  EQ.  3)  THEN 
lOPORT  =  3 

ELSE  IF  (SEL  .  EQ.  4)  THEN 
lOPORT  =  0 

ELSE  IF  (SEL  .  EQ.  5)  THEN 
lOPORT  =50 
ELSE 

lOPORT  =99 
END  IF 

IF  ((SEL  .EQ.  4)  .OR.  (SEL  .  EQ.  5))  THEN 


=>  ') 
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CALL  CLR 
WRITE(*,60) 

WRITE(*,61) 

WRITE(*,62) 

WRITE(*,61) 

WRITE (*,63) 

WRITE(*,61) 

WRITE(*,93) 

60  F0RMAT(lx,///,10x, ' . 

*  . \/,10x,’  ',16x,'BAUD  (data  transfer)  RATE  MENU' 

*,16x,'  ') 

61  FORMATdOx,' . 

*  . d 

62  FORMATdOx,'  '  ,2x, '  SELECT  NO.  '  ,3x, '  '  ,15x, 'BAUD  RATE' ,22x, '  ') 

63  FORMATdOx,'  ',3x,'  1  ',2x,'  ',18x,'  300’ 

*,24x,'  ' 


*,/,10x,' 
*,24x,'  ' 

',3x,' 

2 

’.2x,' 

’ ,18x,' 1200' 

*,/,10x,' 
*,24x,'  ' 

',3x,' 

3 

’,2x,' 

' ,18x,'4800' 

*,/,10x,’ 

'.3x,' 

4 

•,2x,' 

'  ,18x, '9600' 

* , 24x ,  ) 

CALl’cUP( 18,56) 

READ(*,*)  SEL 
IF  (SEL  .EQ.  1)  THEN 
lOPORT  =  lOPORT  -i-  300 
ELSE  IF  (SEL  .EQ.  2)  THEN 
lOPORT  =  lOPORT  +  1200 
ELSE  IF  (SEL  .  EQ.  3)  THEN 
lOPORT  =  lOPORT  +  4800 
ELSE 

lOPORT  =  lOPORT  +  9600 
ENDIF 

CALL  CLR 
WRITE(*,50) 

WRITE(*,51) 

WRITE (*,52) 

WRITE(*,51) 

WRITE(*,53) 

WRITE(*,51) 

WRITE(*,93) 

50  F0RMAT(lx,///,10x, ' . 

*  . ',/,10x,'  '  ,25x,' PARITY  MENU' 

*  26x  '  ' 

51  FORMATdOx,' . 

Vr -  ) 

52  FORMATdOx,'  '  ,2x, '  SE.^'LCT  NO.  '  ,3x, '  ',  17x, 'PARITY' ,23x,  ’  ') 

53  FORMATdOx,'  ',3x,'  1  ',2x,'  ’,14x,'  NO  Parity' 

*,22x,'  ' 

*,/,10x,'  ',3x,'  2  ',2x,’  ' ,14x, 'EVEN  Parity' 

*,2ix,'  ' 

*i/,l6x,'  ',3x,'  3  ',2x,'  ' ,14x, 'ODD  Parity' 

*,22x,'  ') 

CALl’cUP(17,56) 
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READ(*,*)  SEL 
END  IF 

IF  (SEL  .EQ.  3)  THEN 
lOPORT  =  lOPORT  +  1 
ELSE  IF  (SEL  .  EQ.  2)  THEN 
lOPORT  =  lOPORT  +  2 
EENDIF 
99  RETURN 
END 


* 

SUBROUTINE  SIZE  --  ALLOWS  ADJUST  OF  OUTPUT  PLOT  * 

* 


SUBROUTINE  SIZE  (FACT) 

REAL  FACT 

C.  .  .  ADJUST  SIZE  OF  OUTPUT  PLOT 
CALL  CLR 
WRITE (*,15) 

15  F0RMAT(lx,////,10x,' . 

*  - ’,/.10x,  ***  This  is  the  routine  to  Adjust  the  plot  size  ** 

**'./, lOx/ . ' 

*///jlOx,  Output  plots  are  currently  sized  to  fill  the  console  sere 
*en.  ^//, 

*  lOx,  Would  you  like  to  adjust  the  plot  size? ' , lOx, 

*  ' Enter  "y"  or  "n"  ==>  ' ) 

CALL  CUP( 14,31) 

READ(*,'(A)^)  CHG 

IF  ((CHG  .EQ.  'n')  .OR.  (CHG  ,  EQ.  'N'))  THEN 
GO  TO  17 
END  IF 
WRITE(*,16) 

16  FORMAT(/// , lOx, ' A  default  plot  factor  size  of  0.9  is  used  to  fill 
*the  console  screen. lOx, 

*'This  plot  size  will  be  halved  by  entering  a  value  of  0.45. 
*10x,'lt  will  be  doubled  with  a  factor  entry  of  1.8.  ',//,10x, 

*' Enter  factor  size  =>  ') 

CALL  CUP(22,32) 

READ(*,*)  FACT 

17  RETURN 
END 
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APPENDIX  C.  D.S.L  SIMULATION  PROGRAM  AND  DATA  FILE. 


This  appendix  contains  three  listing  files  and  its  resulting  data  files  for  the  D.S.L 
simulation  programs  used  in  examples  4-1,  4-2  and  4-3.  These  programs  are  used  to 
produce  the  data  files  which  generate  the  plots  of  the  performance  index  curves  on  the 
parameter  plane. 


A.  SIMULATION  PROGRAM  FOR  EXAMPLE  4-1. 

TITLE  PERFORMANCE  INDEX  CALCULA"''ON 
ARRAY  Z(1),P(4),C(2),D(2) 

TABLE  Z(l)=100.  ,P(1-4)=1. 0,11. 0, *j.0, 0.0, D(l)=1.0 

PARAM  A=0. 0,B=0. 0,Y0=0.  0 

INITIAL 

C(1)=A 

C(2)=B 

D(2)=B 

BETA=B 

ALPHA=A 

DYNAMIC 

R=STEP(0.  ) 

DERIVATIVE 

Y1=TRNFR(0,3,Y0,Z,P,E) 

M=TRNFR(1,1,Y0,C,D,Y1) 

E=R-M 

ABE=ABS(E) 

PI1=INTGRL(Y0,E*E) 

P I 2=INTGRL( YO , ABE*TIME ) 

PI3=INTGRL(Y0,ABE) 

CONTROL  FINTIM=15,DELT=0.  01 
PRINT  15. , ALPHA, BETA, PI1,PI2,PI3 
END 

PARAM  A=4.  ,B=10. 

END 

PARAM  A=5.  ,B=7.5 
END 

PARAM  A=5.  ,B=10, 

END 

PARAM  A=5.  ,B=15 
END 

PARAM  A=6.  ,B=10 
END 

PARAM  A=7. ,B=30 
END 

PARAM  A=7.  ,B=40 
END 

PARAM  A=7.5,B=15 
END 

PARAM  A=7.5,B=20. 
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END 

PARAM  A=7,5,B=25. 
END 

PARAM  A=7.5,B=35 
END 

PARAM  A=7.5.B=45. 
END 

PARAM  A=7.5,B=50. 
END 

PARAM  A=10.  ,B=10. 
END 

PARAM  A=10.  ,B=15 
END 

PARAM  A=10.  ,B=20. 
END 

PARAM  A=10.  ,B=25. 
END 

PARAM  A=10.  ,B=35 
END 

PARAM  A=10.  ,B=40. 
END 

PARAM  A=12.  5,B=10. 
END 

PARAM  A=12.  5,B=15 
END 

PARAM  A=12.  5,B=20. 
END 

PARAM  A=12.5,B=25. 
END 

PARAM  A=12.5,B=35 
END 

PARAM  A=12.5,B=45. 
END 

PARAM  A=12.5,B=50. 
END 

PARAM  A=15.  ,B=10. 
END 

PARAM  A=15.  ,B=15 
END 

PARAM  A=15,  ,B=20. 
END 

PARAM  A=15.  ,B=25. 
END 

PARAM  A=15. ,B=30. 
END 

PARAM  A=15.  ,B=35 
END 

PARAM  A=15.  ,B=40. 
END 

PARAM  A=15.  ,B=45 
END 

PARAM  A=15. ,B=50. 
END 

PARAM  A=17.5,B=15 
END 

PARAM  A=17.5,B=20. 
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i 

I 


END 


PARAM  A=17.5,B=25. 
END 

PARAM  A=17.5,B=35 
END 


PARAM  A=17.5,B=40. 
END 

PARAM  A=20.  ,B=1S 
END 

PARAM  A=20.  ,B=20. 
END 

PARAM  A=20. ,B=25. 
END 

PARAM  A=20,  ,B=35 
END 

PARAM  A=20. ,B=40. 
END 

PARAM  A=20. ,B=30. 
END 

PARAM  A=20. ,B=60. 
END 

PARAM  A=20. ,B=90 
END 

PARAM  A=25.  ,B=15 
END 

PARAM  A=2S. ,B=20. 
END 


PARAM  A=25. ,B=25. 
END 


PARAM  A=25.  ,B=35 
END 

PARAM  A=25.  ,B=40. 
END 

PARAM  A=25.  ,B=30. 
END 

PARAM  A=30. ,B=20. 
END 

PARAM  A=30.  ,B=25. 
END 

PARAM  A=30. ,B=30 
END 

PARAM  A=30.  ,B=35. 
END 

PARAM  A=30.  ,B=40. 
END 

PARAM  A=50. ,B=20. 
END 

PARAM  A=50. ,B=25. 
END 

PARAM  A=50. ,B=30 
END 

PARAM  A=50. ,B=35. 
END 


PARAM  A=50. ,B=40. 
END 


PARAM  A=50. ,B=50. 
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END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

STOP 


A=50. ,B 
A=60.  ,B 
A=60. ,B 
A=60. ,B 
A=80. ,B 
A=80. ,B 
A=80. ,B 
A=80. ,B 
A=100.  , 
A=100.  , 
A=100.  , 
A=100.  , 
A=100.  , 
A=100.  , 
A=120.  , 
A=120.  , 
A=130.  , 
A=130.  , 
A=140.  , 
A=140.  , 


1=60. 

1=30. 

1=70 

1=90. 

1=30. 

1=50. 

!=70. 

1=90. 

B=30. 

B=40. 

B=50. 

B=60. 

B=80. 

B=100 

B=35. 

B=50. 

B=35. 

B=50. 

B=35. 

B=50. 


71 


B.  DATA  FILE  FOR  EXAMPLE  4-1. 


ALPHA 

BETA 

lES 

lAET 

lAE 

0.0000 

0.000 

1.00000 

1.50000 

1.00000 

4.0000 

10.000 

0.32856 

0.47462 

0.64140 

5.0000 

7.500 

0.30910 

0.53308 

0.64346 

5.0000 

10.000 

0.2S481 

0.33893 

0.54590 

5.0000 

15.000 

0.30S55 

0.42223 

0.60431 

6.0000 

10.000 

0.26177 

0.30259 

0.50802 

7.0000 

30.000 

0.33793 

0.66355 

0.71529 

7.0000 

40.000 

0.41387 

1.1865 

0.92510 

7.5000 

15.000 

0.23340 

0.18931 

0.42445 

7.5000 

20.000 

0.25469 

0.25123 

0.47959 

7.5000 

25.000 

0.28477 

0.38517 

0.56794 

7.5000 

35.000 

0.35289 

0.76717 

0.76028 

7.5000 

45.000 

0.42382 

1.27460 

0.95456 

7.5000 

50.000 

0.45935 

1.57120 

1.05080 

10.000 

10.000 

0.25001 

0.42101 

0.54982 

10.000 

15.000 

0.20401 

0.15825 

0.37735 

10.000 

20.000 

0.21075 

0.14365 

0.37559 

10.000 

25.000 

0.22905 

0.19092 

0.42358 

10.000 

35.000 

0.27603 

0.38136 

0.55835 

10.000 

40.000 

0,30143 

0.50587 

0.62925 

12.500 

10.000 

0.28031 

0.66678 

0.66477 

12.500 

15.000 

0.19186 

0.16546 

0.37150 

12.500 

20.000 

0.18744 

0.11826 

0.33404 

12.500 

25.000 

0.19806 

0.12439 

0.35079 

12.500 

45.000 

0.27154 

0.38180 

0.55436 

12.500 

50.000 

0.29198 

0.47837 

0.61059 

15.000 

10.000 

0.34726 

1.23910 

0.87889 

15.000 

15.000 

0.18889 

0.18998 

0.38491 

15.000 

20.000 

0.17409 

0.11261 

0.31716 

15.000 

25.000 

0.17885 

0.10424 

0.31487 

15.000 

30.000 

0.18996 

0.11504 

0.33643 

15.000 

15.000 

15.000 

15.000 

17.500 

17.500 

17.500 

17.500 

17.500 

20.000 

20.000 

20.000 

20.000 

20.000 

20.000 

20.000 

20.000 

25.000 

25.000 

25.000 

25.000 

25.000 

25.000 

30.000 

30.000 

30.000 

30.000 

30.000 

50.000 

50.000 

50.000 

50.000 

50.000 

50.000 

50.000 


35.000 

40.000 

45.000 

50.000 

15.000 

20.000 

25.000 

35.000 

40.000 

15.000 

20.000 

25.000 

35.000 

40.000 

30.000 

60.000 

90.000 

15.000 

20.000 

25.000 

35.000 

40.000 

30.000 

20.000 

25.000 

30.000 

35.000 

40.000 

20.000 

25.000 

30.000 

35.000 

40.000 

50.000 

6<0(i()0 


0.203S5 

0.21918 

0.23531 

0.25193 

0.19226 

0.16644 

0.16621 

0.1 8434 

0.19669 

0.20112 

0.16250 

0.15766 

0.17018 

0.18022 

0.16203 

0.22763 

0.30424 

0.23759 

0.16216 

0.14800 

0.15137 

0.15796 

0.14718 

0.16949 

0.14448 

0.1 3889 

0.13990 

0.14392 

0.31020 

0.16519 

0.13466 

0.12415 

0. 1 2067 

0.12200 

o.i:’9: 


0.14462 

0.19111 

0.24730 

0.31143 

0.23148 

0.11377 

0.09545 

0.10984 

0.13387 

0.29311 

0.12059 

0.09084 

0.09490 

0.10688 

0.0SSS3 

0.23630 

0.56953 

0.51682 

0.15589 

0.09050 

0.07941 

0.08447 

0.07807 

0.21232 

0.10924 

0.07458 

0.06935 

0.07283 

1.10170 

0.26752 

0.15155 

0.10402 

0.07715 

0.05161 

0.05862 


0.37180 

0.41479 

0.45996 

0.50594 

0.41186 

0.31237 

0.29500 

0.32730 

0.35790 

0.45121 

0.31537 

0.28397 

0.29932 

0.32112 

0.28474 

0.44709 

0.65381 

0.57663 

0.33829 

0.27819 

0.26559 

0.27703 

0.26276 

0.37707 

0.28776 

0.25551 

0.24805 

0.25209 

0.81409 

0.39313 

0.29655 

0.25509 

0.2331  ! 

0,21576 

0.22270 


60.000 

30.000 

0.14256 

0.20934 

0.33390 

60.000 

70.000 

0.12363 

0.05488 

0.21479 

60.000 

90.000 

0.13784 

0.07347 

0.24487 

80.000 

30.000 

0.1 7846 

0.38818 

0.45147 

80.000 

50.000 

0.10841 

0.10385 

0.22883 

80.000 

70.000 

0.10938 

0.05337 

0.19932 

80.000 

90.000 

0.11809 

0.05032 

0.20455 

100.00 

30.000 

0.27014 

0.88639 

0.71319 

100.00 

40.000 

0.12869 

0.23817 

0.31758 

100.00 

50.000 

0.10764 

0.14872 

0.24754 

100.00 

60.000 

0.10224 

0.10517 

0.21883 

100.00 

oO.OOO 

0.10374 

0.06032 

0.19541 

100.00 

100.000 

0.11050 

0.04206 

0.19150 

120.00 

35.000 

0.20936 

0.58018 

0.54931 

120.00 

50.000 

0.11045 

0.19635 

0.26939 

130.00 

35.000 

0.25186 

0.80858 

0.66742 

130.00 

50.000 

0.11297 

0.22142 

0.28185 

140.00 

35.000 

0.32030 

1.26700 

0.85656 

140.00 

50.000 

0.11619 

0.24697 

0.29524 
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r.  SIMULATION  PROGRAM  FOR  EXAMPLE  4-2. 


TITLE  PERFORMANCE  INDEX  CALCULATION 
ARRAY  A(1),B(4) 

TABLE  A(l)=50.  ,B( 1)=1.  0 ,B(4)=0.  0 
PARAM  K1=0. 2,K2=0. 0,y0=0.  0 
INITIAL 

B(2)=3-I-50*K2 

B(3)=2+50*K1 

BETA=K2 

ALPHA=K1 

DYNAMIC 

R=STEP(0.  ) 

DERIVATIVE 

M=TRNFR(0,3,Y0,A,B,E) 

E=R-M 

ABE=ABS(E) 

PI1=INTGRL(Y0,E*E) 

PI 2=INTGRL( YC , ABE*TIME ) 
PI3=INTGRL(Y0,ABE) 

CONTROL  FINTIM=15,DELT=0.  01 
PRINT  15. .ALPHA, BETA, PI1,PI2,PI3 
END 

PARAM  K1=0.  2,K2=0.  1 
END 

PARAM  K1=0. 2,K2=0.  2 
E.ND 

PARAM  K1=0.  2,K2=0.  3 
END 

PARAM  K1=0.  2,K2=0.4 
END 

PARAM  K1=0.  2,K2=0.  5 
END 

PARAM  K1=0.  2,K2=0.  6 
END 

PARAM  K1=0.  2,K2=0.  7 
END 

PARAM  K1=0.  2,K2=0.  8 
END 

PARAM  K1=0.  2,K2=0.  9 
END 

PARAM  K1=0.  2,K2=1.  0 
END 

PARAM  K1=0. 3,K2=0.  1 
END 

PARAM  K1=0. 3,K2=0.  2 
END 

PARAM  K1=0.  3,K2=0.  3 
END 

PARAM  K1=0.  3,K2=0.  4 
END 

PARAM  K1=0. 3,K2=0.  5 
END 

PARAM  K1=0. 3,K2=0.  6 
END 

PARAM  K1=0. 3,K2=0.  7 
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END 

PARAM  K1=0.  3,K2=0.  8 
END 

PARAM  K1=0.  3,K2=0.  9 
END 

PARAM  K1=0.  3,K2=1.  0 
END 

PARAM  K1=0.  4,K2=0.  1 
END 

PARAM  K1=0.4,K2=0. 2 
END 

PARAM  K1=0.  4,K2=0.  3 
END 

PARAM  K1=0.  4,K2=0. 4 
END 

PARAM  K1=0. 4,K2=0. 5 
END 

PARAM  K1=0.4,K2=0.  6 
END 

PARAM  K1=0.  4,K2=0.  7 
END 

PARAM  K1=0.4,K2=0.  8 
END 

PARAM  K1=0.4,K2=0.9 
END 

PARAM  K1=0.  5,K2=0.  1 
END 

PARAM  K1=0.  5,K2=0.  2 
END 

PARAM  K1=0.  5,K2=0.  3 
END 

PARAM  K1=0.5,K2=0.4 
END 

PARAM  K1=0.  5,K2=0.  5 
END 

PARAM  K1=0. 5,K2=0.  6 
END 

PARAM  K1=0. 5,K2=0.  7 
END 

PARAM  K1=0. 5,K2=0.  8 
END 

PARAM  K1=0. 5,K2=0.  9 
END 

PARAM  K1=0.  5,K2=1.  0 
END 

PARAM  K1=0. 6,K2=0. 1 
END 

PARAM  K1=0. 6,K2=0. 2 
END 

PARAM  K1=0. 6,K2=0. 3 
END 

PARAM  K1=0. 6,K2=0. 4 
END 

PARAM  K1=0. 6,K2=0.  5 
END 

PARAM  K1=0. 6,K2=0. 6 
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END 

PARAM  K1=0.  6,K2=0. 7 
END 

PARAM  K1=0.  6,K2=0.  8 
END 

PARAM  K1=0. 6,K2=0.  9 
END 

PARAM  K1=0.  6,K2=1.  0 
END 

PARAM  K1=0.  7,K2=0.  1 
END 

PARAM  K1=0.  7,K2=0.  2 
END 

PARAM  K1=0.  7,K2=0.  3 
END 

PARAM  K1=0.  7,K2=0.4 
END 

PARAM  K1=0.  7,K2=0.  5 
END 

PARAM  K1=0.  7,K2=0.  6 
END 

PARAM  K1=0.  7,K2=0.  7 
END 

PARAM  K1=0.  7,K2=0.  8 
END 

PARAM  K1=0.  7,K2=0.  9 
END 

PARAM  K1=0.  7,K2=1.  0 
END 

PARAM  K1=0.  8,K2=0.  1 
END 

PARAM  K1=0.  8,K2=0.  2 
END 

PARAM  K1=0. 8,K2=0.  3 
END 

PARAM  K1=0. 8,K2=0.  4 
END 

PARAM  K1=0.  8,K2=0.  5 
END 

PARAM  K1=0.  8,K2=0.  6 
END 

PARAM  K1=0. 8,K2=0.  7 
END 

PARAM  K1=0.  8,K2=0.  8 
END 

PARAM  K1=0.  8,K2=0.  9 
END 

PARAM  K1=0.  8,K2=1.  0 
END 

PARAM  K1=0.  9,K2=0.  1 
END 

PARAM  K1=0. 9,K2=0.  2 
END 

PARAM  K1=0. 9,K2=0.  3 
END 

PARAM  K1=0.  9,K2=0.  4 
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END 

PARAM  Kl=i 
END 

PARAM  Kl=i 
END 

PARAM  Kl=( 
END 

PARAM  Kl=l 
END 

PARAM  Kl=i 
END 

PARAM  Kl=( 
END 

PARAM  Kl= 
END 

PARAM  Kl= 
END 

PARAM  Kl= 
END 

PARAM  Kl= 
END 

PARAM  Kl= 
END 

PARAM  Kl= 
END 

PARAM  Kl= 
END 

PARAM  Kl= 
END 

PARAM  Kl= 
END 

PARAM  Kl= 
END 

PARAM  Kl= 
END 

STOP 


. 9,K2=0. 5 
. 9,K2=0.  6 
. 9,K2=0.  7 
. 9,K2=0.  8 
. 9,K2=0.  9 
.  9,K2=1.  0 
. 0,K2=0.  1 
. 0,K2=0.  2 
.  0,K2=0.  3 
. 0,K2=0.  3 
. 0,K2=0.  4 
. 0,K2=0.  5 
. 0,K2=0.  6 
. 0,K2=0.  7 
. 0,K2=0.  8 
.  0,K2=0.  9 
.  6,K2=1.  0 
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D.  DATA  FILE  FOR  EXAMPLE  4-2. 


ALPHA 

BETA 

lES 

lAET 

lAE 

0.20000 

0.00000 

1823.40 

1177.4 

101.49 

0.20000 

0.10000 

0.81565 

4.9245 

1.8653 

0.20000 

0.20000 

0.91714 

6.0067 

2.0813 

0.20000 

0.30000 

1.09550 

8.5450 

2.5026 

0.20000 

0.40000 

1.28850 

11.4640 

2.9349 

0.20000 

0.50000 

1.48580 

14.6030 

3.3547 

0.20000 

0.60000 

1.67890 

17.1180 

3.3763 

0.20000 

0.70000 

1.87450 

20.2920 

4.7094 

0.20000 

0.80000 

2.06100 

22.7820 

4.3853 

0.20000 

0.90000 

2.23130 

24.4710 

4.6228 

0.20000 

1.00000 

2.40750 

26.7370 

4.8862 

0.30000 

0.10000 

0.54210 

1.3228 

1.0492 

0.30000 

0.20000 

0.66416 

2.3429 

1.3567 

0.30000 

0.30000 

0.S02S1 

3.8884 

1.7132 

0.30000 

0.40000 

0.94561 

5.7563 

2.0709 

0.30000 

0.50000 

1.09000 

7.8763 

2.4222 

0.30000 

0.60000 

1.23450 

9.9336 

2.7476 

0.30000 

0.70000 

1.37940 

12.3130 

3.0747 

0.30000 

O.SOOOO 

1.52230 

14.3720 

3.3626 

0.30000 

0.90000 

1.66200 

16.1900 

3.6201 

0.40000 

0.10000 

0.47397 

0.5706 

0.7763 

0.40000 

0.20000 

0.57805 

1.2171 

1.0504 

0.40000 

0.30000 

0.68821 

2.1595 

1.3389 

0.40000 

0.40000 

0.80004 

3.3445 

1.6271 

0.40000 

0.50000 

0.91257 

4.7404 

1.9125 

0.40000 

0.60000 

1.02540 

6.2509 

2.1895 

0.40000 

0.70000 

1.13840 

7.9484 

2.4640 

0.40000 

0.80000 

1.25100 

9.5631 

2.7191 

0.40000 

0.90000 

1.36320 

11.1670 

2.9616 

0.40000 

1.00000 

1.47550 

13.0080 

3.2087 

0.50000 

0.10000 

0.46277 

0.3224 

0.6638 
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0.50000 

0.20000 

0.55073 

0.7599 

0.9005 

0.50000 

0.30000 

0.64156 

1.3845 

1.1390 

0.50000 

0.40000 

0.73322 

2.1813 

1.3762 

0.50000 

0.50000 

0.82524 

3.1372 

1.6120 

0.50000 

0.60000 

0.91744 

4.2318 

1.8455 

0.50000 

0.70000 

1.00970 

5.4501 

2.0760 

0.50000 

0.80000 

1.10200 

6.7009 

2.2984 

0.50000 

0.90000 

1.19420 

8.0218 

2.5160 

0.50000 

1.00000 

1.28650 

9.4955 

2.7345 

0.60000 

0.10000 

0.47534 

0.2533 

0.6421 

0.60000 

0.20000 

0.55087 

0.5490 

0.8301 

0.60000 

0.30000 

0.62798 

0.9905 

1.0287 

0.60000 

0.40000 

0.70557 

1.5579 

1.2275 

0.60000 

0.50000 

0.78336 

2.2462 

1.4262 

0.60000 

0.60000 

0.86125 

3.0514 

1.6245 

0.60000 

0.70000 

0.93921 

3.9528 

1.8211 

0.60000 

0.80000 

1.01720 

4.9258 

2.0147 

0.60000 

0.90000 

1.09520 

5.9872 

2.2067 

0.60000 

1.00000 

1.17320 

7.1377 

2.3975 

0.70000 

0.10000 

0.50008 

0.3878 

0.7400 

0.70000 

0.20000 

0.56605 

0.4620 

0.8138 

0.70000 

0.30000 

0.63298 

0.7805 

0.9736 

0.70000 

0.40000 

0.70021 

1.2014 

1.1408 

0.70000 

0.50000 

0.76757 

1.7177 

1.3100 

0.70000 

0.60000 

0.83499 

2.3256 

1.4801 

0.70000 

0.70000 

0.90245 

3.0156 

1.6500 

0.70000 

O.StxOOO 

0.96993 

3.7809 

1.8192 

0.70000 

0.90000 

1.03740 

4.6278 

1.9882 

0.70000 

1.00000 

1.10490 

5.5291 

2.1552 

O.SOOOO 

0.10000 

0.53189 

0.5459 

0.8400 

0.80000 

0.20000 

0.59036 

0.4757 

0.8496 

O.SOOOO 

0.30000 

0.64946 

0.6791 

0.9597 

0.80000 

0.40000 

0.70875 

0.9953 

1.0967 

O.SOOOO 

0.50000 

0.76813 

1.3928 

1.2406 

0.80000 

0.60000 

0.82756 

1.8649 

1.3873 
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0.80000 

0.70000 

0.80000 

0.80000 

0.80000 

0.90000 

0.80000 

1.00000 

0.90000 

0.10000 

0.90000 

0.20000 

0.90000 

0.30000 

0.90000 

0.40000 

0.90000 

0.50000 

0.90000 

0.60000 

0.90000 

0.70000 

0.90000 

0.80000 

0.90000 

0.90000 

0.90000 

1.00000 

1.0000 

0.10000 

1.0000 

0.20000 

1.0000 

0.30000 

1.0000 

0.30000 

1.0000 

0.40000 

1.0000 

0.50000 

1.0000 

0.60000 

1.0000 

0.70000 

1.0000 

0.80000 

1.0000 

0.9U000 

1.6000 

1.00000 

0.88701 

2.4071 

0.94648 

3.0174 

1.00600 

3.6918 

1.06540 

4.4141 

0.56816 

0.7239 

0.62062 

0.6242 

0.67352 

0.6623 

0.72655 

0.8877 

0.77964 

1.1953 

0.83276 

1.5688 

0.88590 

2.0029 

0.93905 

2.4950 

0.99221 

3.0396 

1.04540 

3.629S 

0.60743 

0.9219 

0.65498 

0.8218 

0.70284 

0.7427 

0.70284 

0.7427 

0.75080 

0.8598 

0.79880 

1.0886 

0.84683 

1.3838 

0.89487 

1.7344 

0.94292 

2.1352 

0.99097 

2.5824 

1.14690 

2.0727 

» 


1.5352 

1.6S38 

1.8324 

1.9801 

0.9400 

0.9400 

0.9825 

1.0862 

1.2059 

1.3317 

1.4606 

1.5913 

1.7226 

1.8540 

1.0400 

1.0400 

1.0448 

1.0448 

1.1062 

1.1999 

1.3055 

1.4169 

1.5314 

1.6476 

1.7147 
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E.  SIMULATION  PROGRAM  FOR  EXAMPLE  4-3. 

TITLE  PERFORMANCE  INDEX  CALCULATION 
ARRAY  A(1),B(3),C(2),D(2) 

TABLE  A(1)=0. 1084,B(l-3)=90. 36,1.  0,0.  0,D(1)=1. 0 
PARAM  Z=0. 0,P=0.  0,Y0=0.  0 
CONST  K1=0. 454616, R=16. 296 
INITIAL 

C(1)=Z 

C(2)=P 

D(2)=P 

BETA=P 

ALPHA=Z 

DYNAMIC 

R=STEP(0.  ) 

DERIVATIVE 

SI=TRNFR(0,2,Y0,A,B,DEL) 

DEL=TRNFR( 1 , 1 , YO , C , D , XI ) 

SIE=SI-R 

X1=K1*SIE 

SIES=SIE*SIE 

DELS=DEL*DEL 

Q=R*SIES+DELS 

J=INTGRL(YO,Q) 

CONTROL  FINTIM=15.  ,DELT=0.  05 
PRINT  15.  , ALPHA, BETA, J 
END 

PARAM  Z=2.  ,P=0.  1 
END 

PARAM  Z=2.  ,P=0.  3 
END 

PARAM  Z=2.  ,P=0.  5 
END 

PARAM  Z=2.  ,P=0. 7 
END 

PARAM  Z=2.  ,P=0.  9 
END 

PARAM  Z=2.  ,P=1.  0 
END 

PARAM  Z=4.  5,P=0.  1 
END 

PARAM  Z=4. 5,P=0.3 
END 

PARAM  Z=4. 5,P=0.  5 
END 

PARAM  Z=4.  5,P=0.  7 
END 

PARAM  Z=4.5,P=0.  9 
END 

PARAM  Z=4.  5,P=1.0 
END 

PARAM  Z=7.  5,P=0.  1 
END 

PARAM  Z=7.5,P=0.3 
END 

PARAM  Z=7.5,P=0.  5 
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END 

PAP  AM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 

END 

PARAM 


Z=7.  5,P=0.  7 
Z=7.  5,P=0.  9 
Z=7,  5,P=1.  0 
Z=9.  ,P=0.  1 
Z=9.  ,P=0.  3 
Z=9. ,P=0. 5 
Z=9. ,P=0. 7 
Z=9. ,P=0. 9 
Z=9.  ,P=1.  0 
Z=13.  ,P=0.  1 
Z=13.  ,P=0.  3 
Z=13.  ,P=0.  5 
Z=13.  ,P=0. 7 
Z=13.  ,P=0.  9 
Z=13.  ,P=1.  0 
Z=16. 5,P=0.  1 
Z=16.  5,P=0.  3 
Z=16. 5,P=0.  5 
Z=16. 5,P=0.  7 
Z=16.  5,P=0.  9 
Z=16. 5,P=1.  0 
Z=19.  ,P=0.  1 
Z=19.  ,P=0.  3 
Z=19.  ,P=0.  5 
Z=19.  ,P=0.  7 
Z=19.  ,P=0.  9 
Z=19.  ,P=1.  0 
Z=23.  5,P=0.  1 
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END 

PARAM  Z=23.5,P=0.  3 
END 

PARAM  Z=23.  5,P=0.5 
END 

PARAM  Z=23.  5,P=0.  7 
END 

PARAM  Z=23.5,P=0.  9 
END 

PARAM  Z=23.5,P=1.0 
END 

PARAM  Z=26.  ,P=0.  1 
END 

PARAM  Z=26.  ,P=0.  3 
END 

PARAM  Z=26.  ,P=0.  5 
END 

PARAM  Z=26.  ,P=0.  7 
END 

PARAM  Z=26.  ,P=0.  9 
END 

PARAM  Z=26.  ,P=1.  0 
END 

PARAM  Z=29.5,P=0.  1 
END 

PARAM  Z=29.  5,P=0.  3 
END 

PARAM  Z=29.5,P=0.5 
END 

PARAM  Z=29.5,P=0.  7 
END 

PARAM  Z=29.5,P=0.  9 
END 

PARAM  Z=29.5,P=1.0 
END 

PARAM  Z=35.  ,P=0.  1 
END 

PARAM  Z=35.  ,P=0. 3 
END 

PARAM  Z=35.  ,P=0. 5 
END 

PARAM  Z=35.  ,P=0. 7 
END 

PARAM  Z=35.  ,P=0. 9 
END 

PARAM  Z=35.  ,P=1. 0 
END 

STOP 
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F.  DATA  FILE  FOR  EXAMPLE  4-3. 


ALPHA 

BETA 

J 

0.0000 

0.00000 

15.000 

2.0000 

0.10000 

14.937 

2.0000 

0.30000 

15.050 

2.0000 

0.50000 

15.042 

2.0000 

0.70000 

15.020 

2.0000 

0.90000 

15.002 

2.0000 

1.00000 

14.994 

4.5000 

0.10000 

14.899 

4,5000 

0.30000 

14.9S2 

4.5000 

0.50000 

14.982 

4.5000 

0.70000 

14.970 

4.5000 

0.90000 

14.959 

4.5000 

1.00000 

14.955 

7.5000 

0.10000 

14.862 

7.5000 

0.30000 

14.912 

7.5000 

0.50000 

14.918 

7.5000 

0.90000 

14.912 

7.5000 

1.00000 

14.911 

9.0000 

0.10000 

14.846 

9.0000 

0.30000 

14.881 

9,0000 

0.50000 

14.889 

9.0000 

0.70000 

14.890 

9.0000 

0.90000 

14.890 

9.0000 

1.00000 

14.890 

13.0000 

0.10000 

14.812 

13.0000 

0.30000 

14.809 

13.0000 

0.50000 

14.820 

13.0000 

0.70000 

14.829 

13.0000 

0.90000 

14.836 

13.0000 

1.00000 

14.839 

16.5000 

0.10000 

14.791 
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16.5000 

0.30000 

14.759 

16.5000 

0.50000 

14.769 

16.5000 

0.70000 

14.7S2 

16.5000 

0.90000 

14.794 

16.5000 

1.00000 

14.799 

19.0000 

0.10000 

14.779 

19.0000 

0.30000 

14.730 

19.000 

0.50000 

14.737 

19.000 

0.70000 

14.752 

19.000 

0.90000 

14.766 

19.000 

1.00000 

14.773 

23.500 

0.10000 

14.765 

23.500 

0.30000 

14.686 

23.500 

0.50000 

14.689 

23.500 

0.70000 

14.704 

23.500 

0.90000 

14.722 

23.500 

1.00000 

14.730 

26.000 

0.10000 

14.760 

26.000 

0.30000 

14.667 

26  000 

0.50000 

14.666 

26.000 

0.70000 

14.681 

26.000 

0.90000 

14.699 

26.000 

1 .00000 

14.708 

29.500 

0. 10000 

14.755 

29.500 

0.300O0 

14.646 

29.500 

0.50000 

14.638 

29.500 

0.70000 

1J.652 

29.500 

0.90000 

14.671 

29.500 

1.00000 

14.680 

35.000 

1.00000 

14.754 

35.000 

0.30000 

14.621 

35.000 

0.50000 

14,603 

35.000 

0.70000 

14.614 

35.000 

0.90000 

14.632 

35.000 

1 .00000 

14.642 
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APPENDIX  D.  USER'S  MANUAL 


When  we  wish  to  plot  the  constant  cost  curves,  the  following  described  procedures 
should  be  performed  as  indicated.  In  this  section,  let  us  consider  the  specific  problem 
in  e.xample  4- 1 . 


I.  PLOTTING  THE  PARA.METER  PLANE  CURVES. 


The  first  step  is  to  plot  the  parameter  plane  curves  which  are  constant  zeta  curves 
and  constant  omega  curves  as  shown  in  Figure  12  in  Chapter  4.  This  curve  permits  us 
to  choose  the  range  of  two  parameters  that  are  to  be  used  in  the  analysis  and  design. 

2.  GENEICXTING  A  DATA  FILE. 


The  second  is  to  obtain  a  data  file  containing  a  ,  ji  and  performance  inde.x  values. 
Simulation  in  D.S.L  is  used  to  evaluate  the  cost  at  a  large  number  of  points  on  the  pa¬ 
rameter  plane.  The  simulation  results  look  as  follows. 

1  DSl/VS .  RElEftSE  1  . HODiriCATIOII  lEVEl  1  .  1985 


0 

0  >>>>  DSl  SIMUlAIlnll  IMPUT  DATA  <<« 

HUE  PERFORHAIICE  INDEX  CAICUIAIIOII 

lADl  EZ(1)=:1 00., P(l-<1)-1. 0,11.0, 10. 0,0.0,D(1)  =  1.D 

FARAM  A=0.n,0=0.0,YO=0.0 

COlURni  FlllIin  =  15,DElI  =  0.01 

PRINI  15.,ALrilA,BElA,IES,lAET,lAE 

END 

ORKR  KIENGN  =  Z189,  XmiNT  =  ?.Z1 ,  AVAIIABIE  SPACE  LEFT  IN  COHHDN/CURVAL/  ■ 

OYX*  RK5  llllENRAriON  HEINOD  USED 

0  RRR  DSL  OUIPUI  IISIING,  GROUP  1 

F  ,  STARTING  RUN  1  Kkr 

PERFORMANCE  INDEX  CALCULATION 


TIME 

ALPHA 

BETA 

lES 

lAET 

lAE 

0. OOOOOEFOO 

0.  OOOOOE'OO 

O.OOOOOEIOO 

0. OOOOOE+00 

O.OOOOOEFOO 

0. OOOOOEFOO 

Q 

15.000 

o.aooooE«oo 

O.OOOOOETOO 

15.000 

112.50 

15.000 

0 

>>>>  DSL 

SIMULATION 

INPUT  DATA  «« 

PARAM 

END 

A=<i.  ,8  =  10. 

0 

KKN  1 

DSL  OUTPUT  LISTING,  GROUP 

1 

+ 

,  STARTING  RUN  2  kkx 

PERFORHANCE  index 

CALCULATION 

TIME 

ALPHA 

DEIA 

1E5 

lAET 

lAE 

O.OOOOOEFOO 

A. 0000 

10.000 

O.DDOOOE4DO 

O.OOOOOEFOO 

O.OOOOOEFOO 

fl 

15.000 

A  .0000 

10.000 

0.32858 

0.A7A62 

0.6A1A0 

0 


W''e  have  to  edit  the  data  file  as  shown  in  Appendi.x  C  to  remove  unnecessary'  parts. 
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3.  LOADING  A  DATA  FILE. 


The  third  step  is  to  download  this  file  from  the  mainframe  computer  to  a  floppy 
diskette.  This  data  now  on  the  lloppy  diskette  can  be  copied  on  the  IBM  microcomputer 
under  a  designated  filename.  In  this  example,  the  filename  is  called  pil.dat.  Now,  we 
are  prepared  to  plot  the  families  of  performance  index  curves  in  the  same  range  of  pa¬ 
rameter  values  on  the  two  parameter  plane. 

4.  PLOTTING  THE  COST  CURVES. 

To  run  the  performance  index  program,  type  in  the  executable  file  name, 
COS  TFUN.  The  IBM  microcomputers  in  the  NFS  Controls  Lab  are  set  up  under  the 
file  management  system  IDIR.  To  execute  the  program  under  the  this  system,  enter  the 
appropriate  subdirectory  by  using  the  up;down  arrow  keys  to  position  the  cursor  adja¬ 
cent  to  the  subdirectory  name,  then  depress  the  EN'l  ER  key.  Use  this  same  procedure 
to  select  COSTFUN.EXE  and  push  ENTER  to  run  the  program.  The  preparation  step 
for  plotting  the  desired  curves  will  appear. 


Enter  proper  values  after  each  prompt  as  follows: 

l.llow  many  constant  PI  curves  #  =  =  =>6 
2. Enter  PI  values  ,lwgt,ldig  f  1  =  =  =  >0 . 27 , 1 , 2 

2. Enter  PI  values  ,lwgt,ldig  #  2===>0.2,i,i 

2. Enter  PI  values  ,lwgt,ldig  #  3  =  =  =  >0 . 1 5 , 1  ,  2 

2. Enter  PI  values  ,lvrgt,ldig  I  4  =  =  =  >0 . 1 2 , 1 , 2 

2. Enter  PI  values  ,lwgt,ldig  #  5  =  =  =  >0 . 1 1 , 1 , 2 

2.  Enter  PI  values  ,lwgt,ldig  #  6===>0. 1,1,1 

3. Enter  number  of  data  file  points  ===>86 
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Next,  adjust  the  plot  size  to  the  same  paratneter  plane  curves. 


This  is  the  routine  to  Adjust  the  plot  size 


Output  plots  are  currently  sized  to  fill  the  console  screen. 

Would  you  like  to  adjust  the  plot  size? 

Enter  ’’y’’  or  "n"  =  =  >  y 

A  default  plot  factor  size  of  0.9  is  used  to  fill  the  console  screen. 
This  plot  size  will  be  halved  by  entering  a  value  of  0.45. 

It  will  be  doubled  with  b  factor  entry  of  1.8. 

Enter'factor  size  ==>  0.7 

The  last  step  is  to  enter  a  data  file  name. 


***  This  is  the  routine  to  LOAD  data  from  a  file  •** 

File  name  should  not  exceed  8  characters 
File  extension  should  not  exceed  3  characters 

What  is  the  file  name  (fn)  and  extension  (ext)? 

Enter  in  form  fn.ext  (e.g.  PARAPLAN. INP)  ===>  pil.dat 


89 


After  displaying  the  graph  on  the  monitor,  the  program  prompts  you  for  a  hardcopy 
of  the  graph. 

WOULD  YOU  LIKE  HARDCOPY  OF  THESE  PLOTS?  (Y  or  N) 


PRINTER  NO. 


1 
2 

3 

4 

5 

6 

7 

8 
9 

++++++♦+++++++++ 
10 
11 
12 
99 


PRINTER/OUTPUT  HENU 

PRINTER 


Epson  FX-80,  All 
Epson  FX-lOO,  All 
Epson  MX-lOO,  All 
Epson  RX-80,  All 
Epson  HX-80  &  IBM  Printer 
HP  7470A  Graphics  Plotter 
HP  7475A  Graphics  Plotter 
HP  758xB  Series  Plotters 
HP  2686A  Laser  Jet 
+++++4+++++++++++t+++++++++++++++++++++++++++++ 
Graphics  Monitor  (default) 

HARDWARE  Interface  Menu 

Input  PLOT88  Values  for  lOPORT  and  MODEL 
RETURN  to  LOAD  Data 


Enter  integer  number  for  selection  =as>  5 

Reenter  the  data  file  name  for  the  hard  copy.  Figure  13  is  a  performance  index 
curve  for  constraint  jBcit.  In  like  manner,  the  other  two  nonlinear  criterion  of  per¬ 
formance  index  in  terms  of /|E|  t  dt  and  J\E\  di  can  also  be  plotted  by  running  the 
same  program.  Figure  Id  and  Figure  15  are  shown  for  desired  cc’vsiraint /I  £  I and 
f\E\  di,  respectively. 
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